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Abstract 

In many extensions of the Standard Model (SM) flavour changing neutral current (FCNC) pro- 
cesses can be mediated by tree- level heavy neutral scalars and/or pseudo-scalars H^{A^). This 
generally introduces new sources of flavour violation and CP violation as well as left-handed 
(LH) and right-handed (RH) scalar (1 =F 75) currents. These new physics (NP) contributions 
imply a pattern of deviations from SM expectations for FCNC processes that depends only 
on the couplings of H^{A^) to fermions and on their masses. In situations in which a single 
or dominates NP contributions stringent correlations between AF = 2 and AF = 1 
observables exist. Anticipating the Flavour Precision Era (FPE) ahead of us we illustrate 
this by searching for allowed oases in the landscape of a given model assuming significantly 
smaller uncertainties in CKM and hadronic parameters than presently available. To this end 
we analyze AF = 2 observables in ^ — ^ and — systems and rare B and K decays 
with charged leptons in the final state including both left-handed and right-handed scalar cou- 
plings of and A^ to quarks in various combinations. We identify a number of correlations 
between various fiavour observables that could test and distinguish these different scenarios. 
The prominent role of the decays Bg^d in these studies is emphasized. Imposing the 

existing flavour constraints, a rich pattern of deviations from the SM expectations in rare Bg ^ 
decays emerges provided Mh < 1 TcV. NP effects in rare K decays are very small. Neutral 
SM Higgs contributions to rare B and K decays turn out to be negligible once the constraints 
from AF = 2 processes are taken into account. Finally, we point out striking differences 
between the correlations found here and in scenarios in which tree-level FCNC are mediated 
by a new neutral gauge boson Z' . 
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1 Introduction 

The recent discovery of a scalar particle with a mass of 126 GeV opened the gate to the 
unexplored world of scalar particles which could be elementary or composite. While 
we will surely learn a lot about the properties of these new objects through collider 
experiments like ATLAS and CMS, also low energy processes, in particular flavour 
violating transitions, will teach us about their nature. In the Standard Model (SM) 
and in many of its extensions there are no fundamental flavour-violating couplings 
of scalar^ to quarks and leptons but such couplings can be generated through loop 
corrections leading in the case of AF = 1 transitions to Higgs-Penguins (HP) and 
in AF = 2 transitions to double Higgs-Penguins (DHP). However, when the masses 
of the scalar particles are significantly lower than the heavy new particles exchanged 
in the loops, the HP and DHP look at the electroweak scale as flavour violating tree 
diagrams. Beyond the SM such diagrams can also be present at the fundamental level, 
an important example being the left-right symmetric models. From the point of view 
of low energy theory there is no distinction between these possibilities as long as the 
vertices involving heavy particles in a Higgs-Penguin cannot be resolved and to first 
approximation what really matters is the mass of the exchanged scalar and its flavour 
violating couplings, either fundamental or generated at one-loop level. While all this 
can be formulated with the help of effective field theories and spurion technology, we 
find it more transparent to study directly tree diagrams with heavy particle exchanges. 

In a recent paper [l] an anatomy of neutral gauge boson (Z' and Z) couplings to 
quark flavour changing neutral currents (FCNC) has been presented. Anticipating 
the Flavour Precision Era (FPE) ahead of us and consequently assuming significantly 
smaller uncertainties in CKM and hadronic parameters than presently available, it was 
possible to find allowed oases in the landscape of new parameters in these models and to 
uncover stringent correlations between AF = 2 and AF = 1 observables characteristic 
for such NP scenarios. 

The goal of the present paper is to perform a similar analysis for scalar neutral parti- 
cles and to investigate whether the patterns of flavour violation in these two different 
NP scenarios (gauge bosons and scalars) can be distinguished through correlations be- 
tween quark flavour observables. Already at this stage it is useful to note the following 
differences in NP contributions to quark flavour observables in these two scenarios: 

• While the lower bounds on masses of Z' gauge bosons from collider experiments 
are at least 1—2 TeV, new neutral scalars with masses as low as a few hundred 
GeV are not excluded. 

• While in the Z' scenarios in addition to new operators also SM operators with 
modified Wilson coefficients can be present, in the case of tree-level scalar ex- 
changes all effective low energy operators are new. 



Unless otherwise specified we will use the name scalar for both scalars and pseudo-scalars. 
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• While there is some overlap between operators contributing to AF = 2 processes 
in Z' and scalar cases after the inclusion of QCD corrections, their Wilson coef- 
ficients are very different. Moreover, in AF = 1 transitions there is no overlap 
with the operators present in Z' models. 

• Concerning flavour violating couplings of Z and the SM Higgs h, in the case of 
the Z boson large NP effects, in particular in rare K decays, are still allowed 
but then its effects in AF = 2 processes turn out to be very small [l]. In the 
Higgs case, the smallness of the Higgs coupling to muons and electrons precludes 
any visible effects from tree-level Higgs exchanges in rare K and B decays with 
muon or electron pair in the final state once constraints from AF = 2 processes 
are taken into account. Simultaneously tree- level Higgs contributions to AF = 2 
transitions can still provide in principle solutions to possible tensions within the 
SM. 

• At first sight the couplings of scalars to neutrinos look totally negligible but if 
the masses of neutrinos are generated by a different mechanism than coupling 
to scalars, like in the case of the see-saw mechanism, it is not a priori obvious 
that such couplings in some NP scenarios could be measurable. Our working 
assumption in the present paper will be that this is not the case. Consequently 
NP effects of scalars in — )■ vr+i/z/, — )• n^uu and b — )• suu transitions will be 
assumed to be negligible in contrast to Z' models, where NP effects in these decays 
could be very important [l]. As we will see, scalar contributions to Kl — )■ 

and Kl — )■ TT^i~^i~ although in principle larger than for K'^ — )■ n^uu, — t- ti^vv 
and h — svv transitions, are found to be small. 

In order to have an easy comparison with the anatomy of FCNCs mediated by a neutral 
gauge bosons presented in [l] the structure of the present work will be similar to the 
structure of the latter paper but not identical, as rare K decays play in this paper 
a subleading role so that emphasis will be put on Bg and Bd systems. In Section [2] 
we describe our strategy by defining the relevant couplings and listing processes to be 
considered. Our analysis will only involve processes which are theoretically clean and 
have simple structure. Here we will also introduce a number of different scenarios for the 
scalar couplings to quarks thereby reducing the number of free parameters. In Section [3] 
we will first present a compendium of formulae relevant for the study of AF = 2 
processes mediated by tree-level neutral scalar exchanges including for the first time 
NLO QCD corrections to these NP contributions. In Section|4]we discuss rare B decays, 
in particular Bg^ — )■ In Section |5] rare K decays are considered. In Section [6] we 

present a general qualitative view on NP contributions to flavour observables stressing 
analytic correlations between AF = 2 and AF = 1 observables. In Section [7] we 
present our strategy for the numerical analysis and in Section [8] we execute our strategy 
for the determination of scalar couplings in the Bg and B^ systems. We discuss several 
scenarios for them and identify stringent correlations between various observables. We 
also investigate what the imposition of the f/(2)^ flavour symmetry on scalar couplings 
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would imply. In Section [9] we present the results for rare K decays, where NP effects are 
found to be small. We also demonstrate that the contributions of the SM Higgs with 
induced flavour violating couplings, even if in principle relevant for AF = 2 transitions, 
are irrelevant for rare K and d decays. A summary of our main results and a brief 



outlook for the future are given in Section 10 



2 Strategy 

2.1 Basic Model Assumptions 

Our paper is dominated by tree-level contributions to FCNC processes mediated by a 
heavy neutral scalar or pseudoscalar. We use a common name, for them unless 
otherwise specified. When a distinction will have to be made, we will either use 
and for scalar and pseudoscalar, respectively or in order to distinguish SM Higgs 
from additional spin particles we will use the familiar 2HDM and MSSM notation: 
iH,A,h). 

Our main goal is to consider the simplest extension of the SM in which the only new 
particle in the low energy effective theory is a single neutral particle with spin and the 
question arises whether this is possible from the point of view of an underlying original 
theory. If the scalar in question is not a SU{2) singlet, then it must be placed in a 
complete SU{2) multiplet, e.g. a second doublet as is the case of 2HDM or the MSSM. 
However, this implies the existence of its SU{2) partners in a given multiplet with 
masses close to the masses of our scalar. In fact in the decoupling regime in 2HDM 
and MSSM the masses of {H^,H^,A^) are approximately degenerate. While SU{2) 
breaking effects in the Higgs potential allow for mass splittings, they must be of 0{v) 
at most and consequently the case of the dominance of a single scalar is rather unlikely. 

It follows then that our scalar should be a SU{2) singlet. In this case, the scalar-quark 
couplings of H come from the following low energy effective operator 

L = Xt^q}iqihsM + h.c. (1) 

with A denoting the cut-off scale of the low energy theory. After the spontaneous 
breakdown of SU{2) the scalar left-handed coupling is given by 

^liH') = ~X1, (2) 
with analogous expression for the right-handed coupling. 

Now, in the case of AF = 2 transitions the scalar contributions are governed only by 
the couplings A^^(iJ°) to quarks and the corresponding Feynman rule has been shown 
in Fig[l] Here (i, j) denote quark flavours. Note the following important property 

A^liH') = [A]j(i/°)]* (3) 
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that distinguishes it from the corresponding gauge couphngs in which there is no chi- 
rahty flip. 

The couphngs A^^(i7°) are dimensionless quantities but as these are scalar and not 
gauge couplings they can involve ratios of quark masses and the electroweak vacuum 
expectation value v or other mass scales. While from the SM, 2HDM and MSSM we are 
used to having scalar couplings proportional to the masses of the participating quarks, 
it should be emphasized that this is not a general property. It applies only if the scalar 
and the SM Higgs, responsible for SU{2) breakdown, are in the same SU{2) multiplet 
or a multiplet of a larger gauge group G. Then after the breakdown of G to SU (2), the 
scalar appears as a singlet of SU (2) symmetry, with couplings to quarks involving their 
masses after SU{2) breakdown. While this is the case in several models, in our simple 
extension of the SM, it is more natural to think that the involved scalar couplings are 
unrelated to the generation of quark masses. 

In spite of the last statement is useful to recall how the quark masses could enter the 
scalar couplings. Which quark masses are involved depends on the model. Considering 
for definiteness the Bg system let us just list a few cases encountered in the literature: 

• In models with MFV in which the scalar couplings are just Yukawa couplings one 
has 

A^^(i/°)oc^, A^^(i7°)oc^ (MFV) (4) 

implying that A^* dominate in these scenarios. Note however that using (jsj) these 
relations also give 

Af(i^°)oc^, At{H')oc'^ (MFV) (5) 

which implies some care when stating whether LH or RH scalar couplings are 
dominant. Below we will use the ordering bq for AF = 2 operators in Bq (q = s, d) 
systems while qb in the case of rare Bq decays. In the K system sd will be used 
for both AF = 2 and AF = 1 couplings. 

• In non-MFV scenarios the mass dependence in scalar couplings can be reversed 

A5^^(/70) oc ^, A^^(i7°)oc^ (non-MFV) (6) 



implying that A^ dominate in these scenarios. Correspondingly ([5j) is changed to 
Af(i7°)oc^, A?j^(i7°)oc^ (non-MFV) (7) 

promoting the so-called primed operators in AF = 1 decays. 

There exist also models in which flavour violating neutral scalar couplings do not 
involve the masses of external quarks. This is the case for the neutral heavy Higgs 
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in the left-right symmetric models analysed in [2] where the scalar down quark 
couplings are proportional to up-quark masses, in particular rrtt/v. In the case of 
a manifest left-right symmetry with the right-handed mixing matrix being equal 
to the CKM matrix one finds 



equal to the CKM matrix is ruled out by the data, there could be other model 
constructions in which ^ could be satisfied. Also the LH and RH couplings 
differing by sign could in principle be possible. 



Figure 1: Feynman rules for a neutral colourless scalar particle with mass Mh, 
where i, j denote different quark flavours and a, /3 the colours. Pl,r = (1 T 75)/2- 



2.2 Scenarios for Scalar Couplings 

In order to take these different possibilities into account and having also in mind that 
scalar couplings could be independent of quark masses, we consider the following four 
scenarios for their couplings to quarks keeping the pair {i,j) fixed: 

1. Left-handed Scenario (LHS) with complex A^^ and A^^ = 0, 

2. Right-handed Scenario (RHS) with complex A^ ^ and A^^ = 0, 

3. Left-Right symmetric Scenario (LRS) with complex A^^ = Ag ^ 0, 

4. Left-Right asymmetric Scenario (ALRS) with complex A^^ = -Ag ^ 0, 

with analogous scenarios for the pair {s,d). For rare Bg decays in which the ordering 
qb is used, the rule ^ has to be applied to each scenario. For K physics this is 
not required. In the course of our paper we will list specific examples of models that 
share the properties of these different scenarios. We will see that these simple cases 
will give us a profound insight into the flavour structure of models in which NP is 
dominated by left-handed scalar currents or right-handed scalar currents or left-handed 
and right-handed scalar currents of the same size. We will also consider a model in which 
both a scalar and a pseudoscalar with approximately the same mass couple equally to 
quarks and leptons. Moreover we will study a scenario with underlying flavour f/(2)^ 
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symmetry which will imply relations between A^f and A^f couplings and interesting 
phenomenological consequences. 

The idea of looking at NP scenarios with the dominance of certain quark couplings to 
neutral gauge bosons or neutral scalars is not new and has been motivated by detailed 
studies in concrete models like supersymmetric flavour models |3|, LHT model with 
T-parity j4|[5] or Randall-Sundrum scenario with custodial protection (RSc) [6j. See 
also (7||8]. Also our recent analysis of tree- level FCNCs mediated by Z' and Z in [l] 
demonstrates this type of NP in a transparent manner. 

2.3 Scalar vs Pseudoscalar 

It will turn out to be useful to exhibit the differences between the scalar and pseu- 
doscalar spin particles, although one should emphasize that in the presence of CP 
violation, the mass eigenstate propagating in a tree-diagram is not necessarily a CP 
eigenstate. Therefore, generally the coupling to [i^ ii~ appearing at many places in our 
paper can have the general structure 

L = \mf{ii')^irAfm)ii (9) 

where generalizing the Feynman rule in Fig. [T] to charged lepton couplings we have 
introduced: 

Af (i/) = + Af (i/), 

(10) 

t^f[H) = A'^J'iH) - Af{H). 

Ag'^ is real and Ap'^ purely imaginary as required by the hermiticity of the Hamiltonian 
which can be verified by means of ([3]). 

The expressions for various observables will be first given in terms of the couplings 
A^l j^{H) and A'g^p^H) and can be directly used in the case of the scalar particle being 
CP-even eigenstate, like {H°, h) in 2HDM or MSSM setting Af{H) = 0. However, 
when the mass eigenstate is a pseudoscalar A, implying A^^^^ = 0, it will be useful to 
exhibit the i which we illustrate here for the system: 

Al^{A) = -^Ai^{A), A'iiA) = +zA??(A), Af{A) = ^Af{A). (11) 

Here the flavour violating couplings A^f j:j(A) are still complex, while A'^{A) is real. 
The following useful relations follow from ([S]) and (11): 



A^^iA) = z[A'j^{A)]\ Af{A) = -^[A''i{A)r. (12) 



As far as AF = 2 transitions are concerned this distinction between scalar and pseu- 
doscalar mass eigenstate is only relevant in a concrete model in which the relevant 



2 Strategy 



7 



couplings are given in terms of fundamental parameters. However, as in our numerical 
analysis we will treat the flavour violating quark-scalar couplings as arbitrary complex 
numbers to be bounded by AF = 2 observables it will not be possible to distinguish a 
scalar and pseudoscalar boson on the basis of AF = 2 transitions alone. On the other 
hand, when rare decays, in particular Bs^d — ^ i are considered there is a difference 
between these two CcLSGS cLS the pseudoscalar contributions interfere with SM contribu- 
tion, while the scalar ones do not. Consequently the allowed values for /S.^p{A) and 
A'il^{A) will differ from each other and we will find other differences. Finally, if both 
scalar and pseudoscalar contribute to tree-level decays and have approximately the same 
mass as well as couplings related by symmetries, also their contributions to AF = 2 
processes differ. We will consider a simple example in the course of our presentation. 



2.4 Steps 

Let us then outline our strategy for the determination of flavour violating H couplings 
to quarks and for finding correlations between flavour observables in the context of the 
simple scenarios listed above. Our strategy will only be fully effective in the second half 
of this decade, when hadronic uncertainties will be reduced and the data on various 
observables significantly improved. It involves ten steps including a number of working 
assumptions: 

Step 1: 

Determination of CKM parameters by means of tree-level decays and of the necessary 
non-perturbative parameters by means of lattice calculations. This step will provide the 
results for all observables considered below within the SM as well as all non-perturbative 
parameters entering the NP contributions. As \Vub\ is presently poorly known, it will 
be interesting in the spirit of our recent papers [l|[2|[9] to investigate how the outcome 
of this step depends on the value of \Vub\ with direct implications for the necessary size 
of NP contributions which will be different in different observables. 

Step 2: 

We will assume that the ratios 

for scalar and pseudoscalar bosons have been determined in pure leptonic processes and 
that the scalar couplings to neutrinos are negligible. The properties of these couplings 
have been discussed above. In principle these ratios can be determined up to the sign 
from quark flavour violating processes and in fact we will be able to bound them from the 
present data on Bs — )■ but their independent knowledge increases predictive power 

of our analysis. In particular the knowledge of their signs allows us to remove certain 
discrete ambiguities and is crucial for the distinction between LHS and RHS scenarios 
in fJ'^f^^^ decays. Of course, in concrete models like 2HDM or supersymmetric 

models these couplings depend on the fundamental parameters of a given model. 
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Step 3: 

Here we will consider the system and the observables 

where A'^jt^ and can be extracted from the time-dependent Bs — > rate 
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11 . Explicit expressions for these observables in terms of the relevant couplings can oe 



found in Sections [3] and HI 

Concentrating in this step on the LHS scenario, NP contributions to these three ob- 
servables are fully described by 

g^3-...3 ^rp(^) .5. 

with the second ratio known from Step 2. Here §23 > and it is found to be below 
unity but it does not represent any mixing parameter as in 12 . The minus sign is 
introduced to cancel the minus sign in Vts in the phenomenological formulae listed in 
the next section. 

Thus we have five observables to our disposal and two parameters in the quark sector to 
determine. This allows to remove certain discrete ambiguities, determine all parameters 
uniquely for a given Mh and predict correlations between these five observables that 
are characteristic for this scenario. 

Step 4: 

Repeating this exercise in the B*^ system we have to our disposal 

AM,, S^Ks. i3(i?,^/iV), ^'v- (16) 

Explicit expressions for these observables in terms of the relevant couplings can be 
found in Sections [3] and HI 

Now NP contributions to these three observables are fully described by 

Mh Mh ' Mh ' 

with the last one known from Step 2 and bounded in Step 3. Again we can determine 
all the couplings uniquely for a given Mh- Our notations and sign conventions are as 
in Step 3 with S13 > but no minus sign as Vtd has no such sign. 

Step 5: 

Moving to the K system we have to our disposal 

Sk, KL-^Ti^'tt, KL^fi^fi-, (18) 
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where in view of hadronic uncertainties the last decay on this hst will only be used 
to make sure that the existing rough bound on its short distance branching ratio is 
satisfied. Unfortunately tree- level neutral Higgs contributions to — )■ tt^vv and 
Kl — TT^vv are expected to be negligible, but this fact by itself offers an important 
test and distinction from tree level neutral gauge boson exchanges where these decays 
could still be significantly affected (Tj. Also the decays Kl — )• ti^I~^1~ are subject to 
considerable hadronic uncertainties and their measurements are not expected in this 
decade. Yet, as they are known to be sensitive to NP effects it is of interest to consider 
them as well and compare the scalar case with the case of Z' models \\\. 

In the present paper we do not study the ratio e' /e, which is rather accurately measured 
but presently subject to much larger hadronic uncertainties than observables listed in 



(18). Yet, it should be emphasized that e' /e is important for the tests of H FCNC 



scenarios as it is very sensitive to any NP contribution 13 -15 . 

Explicit expressions for the observables in the K system in terms of the relevant cou- 
plings can be found in Sections [3] and [5j 

Now NP contributions to these observables are fully described by 

Mh Mh ' Mh 

The ratios involving muon couplings are already constrained or determined in previous 
steps. Consequently, we can bound quark couplings involved by using the data on the 



observables in (18). Moreover we identify certain correlations characteristic for LHS 
scenario. Si2 > and the minus sign is chosen to cancel the one of Vts- 

We can already announce at this stage that the results in K physics turned out to be 
much less interesting than in the Bg and systems and we summarize them separately 
in Section [9l 

Step 6: 

As all parameters of LHS scenario have been fixed in the first five steps we are in the 
position to make predictions for the following processes 

B^xj+r, B^Ki+r, B^K*tr (20) 

and test whether they provide additional constraints on the couplings. Again as in 
the case of — )■ vr+i/z/ and Kl — )■ 'k^vv also the h — ?• svv transitions are expected 
to be SM-like which provides a distinction from the gauge boson mediated tree-level 
transitions [l]. 

Step 7: 

We repeat Steps 3-6 for the case of RHS. We will see that in view of the change of 
the sign of NP contributions to Bg ^ — )• and — )• ii^ji" decays the structure of 

the correlations between various observables will distinguish this scenario from the LHS 
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one. Yet, as we will find out, by going from LHS to RHS scenario we can keep results 
of Steps 3-5 unchanged by interchanging simultaneously two big oases in the parameter 
space that we encountered already in our study of the 331 model 12 and Z' models [1]. 



This LH-RH invariance present in Steps 3-5 can be broken by the b — >■ si~^i transition 



in (20). They allow us to distinguish the physics of RH scalar currents from LH ones. 



As only RH couplings are present in the NP contributions in this scenario, we can use 



the parametrization of these couplings as in (15), (17) and (19) keeping in mind that 
now RH couplings are involved. 

Step 8: 

We repeat Steps 3-6 for the case of LRS. In the case of tree-level gauge boson contribu- 
tions the new features relative to the previous scenarios is enhanced NP contributions 
due to the presence of LR operators in AF = 2 transitions. Yet, in the scalar case, 
the matrix elements of SLL and SRR operators present in previous scenarios are also 
significant larger than the SM ones and the addition of LR operators has a more modest 
effect than in the gauge boson case. However, one of the important new feature is the 
vanishing of NP contributions to Bg^d — ^ fJ'~^f^~ and Kl — ?• fJ'~^fJ'~ decays. As the LH 
and RH couplings are equal we can again use the parametrization of these couplings 



as in (15), (17) and (19) but their values will change due to different constraints from 
AF = 2 transitions. Also in this step b — )■ si'^i" transitions can play an important 
role. 

Step 9: 

We repeat Steps 3-6 for the case of ALRS. Here the new feature relatively to LRS are 



non-vanishing NP contributions to Bs^d , including 5"^+ _ CP asjTumetries. 

Again the b — > si'^i' transitions will exhibit their strength in testing the theory in a 
different environment: NP contributions to AF = 2 observables due to the presence of 
LR operators. As the LH and RH couplings differ only by a sign we can again use the 



parametrization of these couplings as in (15), (17) and (19) but their values will change 
due to different constraints from AF = 2 transitions. 

Step 10: 

One can consider next the case of simultaneous LH and RH couplings that are unrelated 
to each other. This step is more challenging as one has more free parameters and in 
order to reach clear cut conclusions one would need a concrete model for H couplings 
or a very involved numerical analysis iTjlS, 16 . A simple model in which both a scalar 



and a pseudoscalar with approximately the same mass couple equally to quarks and 



leptons has been recently presented in 17 showing that the structure of correlations 



can be quite rich. We refer to this paper for details. 

Once this analysis of H contributions is completed it will be straightforward to apply 
it to the case of the SM Higgs boson with flavour violating couplings. Yet, we will see 
that this case is less interesting than the case of Z with flavour violating couplings. 
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3 AF = 2 Processes 



3.1 Preliminaries 

In the SM the dominant top quark contributions to AF = 2 processes are described by 
flavour universal real valued function given as follows (xj = m^/M^): 

ixt - Ux^t + xl ?,xl log Xt , . 

- 4(1 -a:,)^ - 2(13^ " ^'^^ 

In other CMFV models SqIxi) is replaced by a different function which is still flavour 
universal and is real valued. This implies very stringent relations between various 
observables in three meson system in question which have been reviewed in (Tsl. 



In the presence of H tree-level contributions the flavour universality is generally broken 
and one needs three different functions 

S{K), S{B,), S{Bs), (22) 

to describe — and B^^ — B^^ systems. Moreover, they all become complex quan- 
tities. Therefore CMFV relations are generally broken. In introducing these functions 
we will include in their definitions the contributions of operators with LL, RR and LR 
Dirac structures. 

The derivation of the formulae listed below is so simple that we will not present it here. 
In any case, the compendium of relevant formulae given below and in next sections is 
self-contained as far as numerical analysis is concerned. 



3.2 Master Functions Including H Contributions 

Calculating the contributions of H to AF = 2 transitions it is straightforward to write 



down the expressions for the master functions S{M) in (22) in terms of the couplings 
defined in Fig. [T} 

We define first the relevant CKM factors 

Xf^ = V*V,,, Xf^ = V,lVu, X['^ = V:,Vts, (23) 

and introduce 

91m = ^2n = 1-78137 X IQ-^GeV-^ . (24) 

a/2 2tt sm^ Uw 

The AF = 2 master functions for M = K, Bg are then given as follows 

S{M) = Soixt) + A5(M) = \S{M)\e''^s' (25) 
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with AS{M) receiving contributions from various operators so that it is useful to write 
AS{M) = [A5(A'/)]sLL + [A5(M)]sRR + [A5(M)]lr. (26) 



The contributing new operators are defined for the K system as follows 19 ,20] 

Qj^^ = {s^.PLd) {sYPRd) , (27a) 



Q^^ = (sPLd) (sPRd) . (27b) 

gf L = (sPLd) (sPLd) , (28a) 

Qf"^ = (sPnd) (rsPnd) , (28b) 

Qf"- = (sa^M) {sa^'PLd) , (28c) 

Qf-"" = {sa,,PRd) (m^-'PRd) , (28d) 

where Pr^l = (1 ± 75)/2 and we suppressed colour indices as they are summed up in 
each factor. For instance s'-^^Pid stands for Sa'jfj.Pida and similarly for other factors. 
For 5° — 5° mixing our conventions for new operators are: 

(29a) 
(29b) 

(30a) 
(30b) 
(30c) 
(30d) 

In order to calculate the SLL, SRR and LR contributions to AS{M) we introduce 
quantities familiar from SM expressions for mixing amplitudes 

^^^'^^ = ^2FBMmBM'^ [xi'^r^B, (31) 
TiK) = ^.FlBKruKM'^ (a^ ^)%2, (32) 





[h.PLq) {prPRq) , 




{hPLq) (bPRq) , 


gf^ = 


(bPLq) {hPLq) , 


gSRR^ 


{hPRq) {bPRq) , 




[ha^^PLq] {ba^'PLq) 


Qf^ = 


{bcr,,PRq) [ha^^'PRq) 
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where r]i are QCD corrections and Bi known SM non-perturbative factors. 
Then 

(33) 

with the SRR contribution obtained by replacing L by R. Note that this replacement 
only affects the coupling Af^{H) as the hadronic matrix elements being evaluated in 
QCD remain unchanged and the Wilson coefficients have been so defined that they also 
remain unchanged. For LR contributions we find 

T{K)[AS{K)U = [c^R^f,^mK^fiH,K))+C^\f,H){Q^^{f^H,K))] 

(34) 



Including NLO QCD corrections 20 the Wilson coefficients of the involved operators 
are given by 

^(/.) = (/.) = 1 + £ (-3 log ^ + , (35) 

Cl--i,) = Cf-(^) = ^ (-^ log ^ + , (36) 

C^r(/^) = -^^, (37) 
CLR(^) = 1-^1 = 1-^. (38) 

Next 

{Q1{f^H,K))^'^^Pt{fiH,K) (39) 
are the matrix elements of operators evaluated at the matching scale fin = 0{Mh) and 



are the coefficients introduced in 19 . The hh dependence of C^djH) cancels the 



one of P°-{ijlh) so that S{K) does not depend on hh- 
Similarly for Bq systems we have 

T(S,)[A5(S,)]sLL = -^^1^ [CfL(/.^)(Qf^(/.H,S,)) ^(/i/.)(Qf"(/^H,5,))] 

(40) 

nB,)[AS{B,)U = -M^Il® [cl^^^^){Q\^(^^^,B,)) + C^^M{Q^^ifiH,B,)^ , 

(41) 
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where the Wilson coefficients C^{^h) are as in the K system and the matrix elements 
are given by 



{Qt{^^H,B,)) 



(42) 



For SRR contributions one proceeds as in the K system. 

Finally, we collect in Table hj central values of {QKhh))- They are given in the MS-NDR 



scheme and are based on lattice calculations in 



21 



22 



for — system and in 23 



for . — o systems. For the — system we have just used the average of the 



results in 21 , 22 that are consistent with each other. As the values of the relevant Bi 
parameters in these papers have been evaluated at /i = 3 GeV and 4.2 GeV, respectively, 
we have used the formulae in |jl9] to obtain the values of the matrix elements in question 
at hh- For simplicity we choose this scale to be Mh but any scale of this order would 
give the same results for the physical quantities up to NNLO QCD corrections that are 
negligible at these high scales. The renormalization scheme dependence of the matrix 
elements is canceled by the one of the Wilson coefficients. 

In the case of tree-level SM Higgs exchanges we evaluate the matrix elements at mt{mt) 
as the inclusion of NLO QCD corrections allows us to choose any scale of 0{Mh) 
without changing physical results. Then in the formulae above one should replace Mh 
by the SM Higgs mass and /ih by mt{mt). The values of hadronic matrix elements at 
mt{mt) in the MS-NDR scheme are given in Table ^ 











{QfH^^H)) 




-0.14 


0.22 


-0.074 


-0.128 


B',-B', 


-0.25 


0.34 


-0.11 


-0.22 


B's-m 


-0.37 


0.51 


-0.17 


-0.33 



Table 1: Hadronic matrix elements {Q^{^h)) in units of GeV^ at fiH = 1 TeV. 















-0.11 


0.18 


-0.064 


-0.107 


B^B^, 


-0.21 


0.27 


-0.095 


-0.191 


B's-B', 


-0.30 


0.40 


-0.14 


-0.29 



Table 2: Hadronic matrix elements (Qiifit)) in units of GeV^ at mt{mt). 
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3.3 Basic Formulae for AF = 2 Observables 

The AB = 2 mass differences are given as follows: 



AM, = ^,M^mB,\XiyFlBB,VB\S{B,)\ , (43) 

AM, = ^,M'^mB^xP\'FlBB^r^B\SiBs)\ • (44) 
The corresponding mixing induced CP-asymmetries are then given by 

S'^i^g = sin(2/3 + 2v9bJ , S^^ = sm{2\(3s\ - 2(fBs) , (45) 
where the phases /3 and (3s are defined by 

Vtci=\Vt,\e-'^, = -\Vts\e-'^\ (46) 

(3s — — 1° . The new phases ipB^ are directly related to the phases of the functions 

s{B,y. 

2^B, = -9 s'- (47) 

Our phase conventions are as in |lj and our previous papers quoted in this work. Con- 
sequently ~ 0.04. On the other hand the experimental results are usually given 
for the phase 

0, = 2Ps + 0'''' (48) 

so that 

5v-<A = -sin(0,), 2(^5, =0NP. (49) 
Using this dictionary the most recent result for 0, from the LHCb analysis of CP- 



violation in Bs — )■ tpfj) decay implies 24 



2|/3,| -2y?B^ = 0.001 ±0.104, (50) 

that is close to its SM value. But the uncertainties are still sufficiently large so that it 
is of interest to investigate correlations of 5*^0 with other observables in the Bs system. 

For the CP-violating parameter ek and AMi^ we have respectively 
where 

(Mi^)* = ^F'^EktukM^ [X%x, + \lmS{K) + 2K\t7^^S,{x,, Xt)] . (52) 
Here, 5*0 ( real valued one-loop box function for which explicit expression is 



given e. g. in 25 . The factors are QCD corrections evaluated at the NLO level 
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in p6}]30] . For rji and rjs also NNLO corrections are known 31 32]. Next 
0.05)° and k _ 



0.94 ± 0.02 

distance effects in 53(ri2 



takes into account that (p^ ^ 



and '^[M_ 



(43.51 ± 
I and includes long 



In the rest of the paper, unless otherwise stated, we will assume that all four parameters 
in the CKM matrix have been determined through tree-level decays without any NP 
pollution and pollution from QCD-penguin diagrams so that their values can be used 
universally in all NP models considered by us. 



4 Rare B Decays 
4.1 Preliminaries 

These decays played already for many years a significant role in constraining NP mod- 
els. In particular Bs — )■ fJ''^fi~ was instrumental in bounding scalar contributions in the 
framework of supersymmetric models and two Higgs doublet models (2HDM). Recently 
a very detailed analysis of the decay Bg — > including the observables involved 



in the time-dependent rate has been presented 17 . Below, after recalling the rele- 
vant effective Hamiltonian that can be used for other b si'^i~ transitions, we will 
summarize the final formulae for the most important observables in Bg — )■ fi~^ii~ that 



have been derived and discussed in more detail in 17 and in particular earlier in 11 
While our analysis of Bs — )■ is less detailed than the one in 17 , our main goal 

here is to discuss the correlations of Bg observables with AF = 2 observables, 

in particular S^^, which were not presented there. Moreover, we analyze here similar 
correlations involving B^ fi~^fi~ observables and S^Ks- 

4.2 Effective Hamiltonian for b — )- si'^i~ 

For our discussion of B^i s — )■ fi^fi^ and for the imposition of the constraints from other 
b si+i- transitions, like B -> K*£+£-, B K£+£- and B Xs£+£-, we will need 
the corresponding effective Hamiltonian which is a generalization of the SM one: 

nes{b^s££)=neE{b^s^)-^^V*V,, J2 [a(/i)g.(/i)+q(^)g:(/i)] (53) 

^ i=9,lO,S,P 

where 





= {sj,PLb){ir£), 


Q'9 = 


{si,PRb){hn), 


(54a) 


Qio 


= {si.PMlrid). 


Qio ~ 


- {-si,PRb){£Yi,£), 


(54b) 


Qs 


= mt{sPnb){££), 


Q's = 


m,{sPLb){££), 


(54c) 


Qp 


= m,{sPRb){h,£), 


Q'p = 


m,{sPLb){H£). 


(54d) 
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Including the factors of nib into the definition of scalar operators makes their matrix 
elements and their Wilson coefficients scale independent. Ti es{b — 57) stands for the 
effective Hamiltonian for the b ^ sj transition that involves the dipole operators. We 
will not discuss 6 — j- S7 in this paper as it appears first at one-loop level and a neutral 
scalar contribution would only be of relevance in the presence of fiavor- conserving scalar 
couplings to down-type quarks which would introduce new parameters without any 
impact on our results. 

Note the difference of ordering of flavours relatively to AF = 2 as already stressed in 
Section |2} Therefore the unprimed operators Qs and Qp represent the LHS scenario 
and the primed ones Q'^ and Q'p the RHS scenario. We neglect effects proportional to 
rris in each case but keep and rrifi different from zero when they are shown explicitly. 

The Wilson coefficients Cg and Cio do not receive any new contributions from scalar 
exchanges and take SM values 

sin^ ewCl"^ = [VYYoixt] - 4sin2 OwZoixt)], (55) 

sm'ewC!^ = -VYYo{xt). (56) 

On the other hand with <^ we have Cg = C[q = 0. Here Yq^Xi) and Zq^Xi) are 
SM one-loop functions given by 

N Xt [xf — A SxjlogxA , „, 

, , 1 ISx^ - 163x3 + 259x2 - 108x 32x^ - SSx^ - 15^^ + 18x , 
Z,{x) = -- logx + 144(^ _ 1)3 + n^^^ • 

(58) 

The coefficient vjY is a QCD factor which for mt = mt{mt) is close to unity: riy = 1.012 



For the coefficients of scalar operators we find 

m^sm c'vi/Os — , [oy) 



^w^s -^:2~YJ2 ^T^TT ' y'^^) 

(61) 



1 


1 








VtlVtb 


1 


1 


AtiH)AfiH) 




M% 


VtlVtb 


1 


1 


A^^{H)Af{H) 




M% 


VtlVtb 


1 


1 


Af{H)Af{H) 



(62) 
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where A'^'^p[H) are defined in (10). We recall that in terms of the couplings used in the 
analysis of S° ^ — 5° ^ mixings we have 

AtiH) = [Al%H)r, At{H) = [A'im*, (63) 

which should be kept in mind when studying correlations between AF = 1 and AF = 2 
transitions. These relations can be directly used in the case of Cs and Cg but in the 
case of Cp and C'p, as discussed in Section [2| it is useful to use in this context the 
following relations: 

A^^^Af = -[A^^]*Af , Af A^^ = [A^^]*Af (64) 
with Ap^ being imaginary but Ap^ real. 



4.3 Observables for 



In the general analysis of Bg — t- fi^fi in 17 , which goes beyond the NP scenario 



considered here, the basic four observables are 



R, 



S\ 



(65) 



Here, the observable i?, defined in (68), is just the ratio of the branching ratio that 



includes AF^ effects and of the SM prediction for the branching ratio that also includes 
them. The relation oi Rio R introduced in 



11 



is given below. Following 17 we will 



denote branching ratios containing AF^ effects with a har while those without these 
effects without it. 

The next two observables, A^^^ '^/^m ^^"^ be extracted from flavour untagged and 
tagged time-dependent measurements of Bs — )• fi'^fi', respectively. As these three 
observables depend also on the new phase (pBs the B^ — B^ mixing, also the mixing 
induced CP-asymmetry S^^ is involved here. 

In order to calculate these observables one introduces 



•-"10 



+ 



Cp — C'p 



2m^ nib + ms 



'-^10 



(66) 



S 



Ami m? 



Bs 



rrib 



Cq — C'q 



m| 2m^ nib + Cf^ 



\S\e 



(67) 



One finds then three basic formulae 111 17l37 



R 



ISM 



1 -|- v4. Us 



l + Vs cos{2ipp - 2ipbJ 



l + Vs 



\pr + 



xi\P\' + \S\') 
l-Vs cos(2v9g - 2v?B. 



(6J 
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OS 



Here 38 







|P|2 + 








|P|2+ |^|2 








l/, = rB,^ = 0.088 ±0.014. 





(69) 



(70) 



(71) 



The ratio R of 



11 



which did not include AF^ effects in the SM result and R which 
includes them are related by 

R = il-ys)R. (72) 

The advantage of R over R is that in the SM it is equal to unity and its departure from 
unity summarizes total NP effects present both in the decay and mixing. 

Another useful variable encountered in this discussion is 



1 - y's 



1 + vA^p'" Us 



(73) 



It is the correction factor that one has to introduce in any model in order to compare 
the branching ratio calculated in this model without AF^ effects and the branching 
ratio which includes them flO 



11 



39 



(74) 



It should be emphasized that presently only B{Bs — )■ /U+yU^) is known experimentally 
but once A^^y ^^^^ extracted from time-dependent measurements, we will be able 
to obtain B{Bs — )■ /x"*"//") directly from experiment as well. Evidently, in any model 
the branching ratios without AFg effect are related to the corresponding SM branching 
ratio through 

B{B, ^ /iV) = B{B, ^ /iV)sM(|P|' + l^n. (75) 

As Fb^ cancels out in the evaluation of A^^y- ^^"^ '^ma" these are theoretically clean 
observables and offer new ways to test NP models. Indeed, as seen in (69) and (70), 
both observables depend on NP contributions and this is also the case of the conversion 
factor r{ys). In the SM and CMFV models P = 1 and S" = so that 

A^J^Y = 1, S'^^ = 0, r{ys) = 0.912 ± 0.014 (SM, CMFV), (76) 

independently of NP parameters present in the whole class of CMFV models. 



As A^Y does not rely on flavour tagging, which is difficult for a rare decay, it will be 
easier to determine than S^^. Given limited statistics, experiments may first measure 
the Bg — )• effective lifetime, a single exponential fit to the untagged rate, from 

which A'^Y ^Iso be deduced [ll]. See also jlT] for discussion. 

While AFrf is very small and yd can be set to zero, in the case of B^ — )■ /^^/^^ one can 



still consider the CP asymmetry 5**^ 



37 



for which one can use all expressions given 



above with the flavour index "s" replaced by "ci" . 
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4.4 Present Data 



The most recent results from LHCb read 40 ,41 



B{Bs fi+fi-) = (3.21};^) X 10"^ B{Bs -> ;U+/i-)sM = (3.25 ± 0.17) x 10-^ (77) 

i3(fid ^/iV") < 9.4 X 10^1°, i3(5d ^ /i+/i")sM = (1.05 ± 0.07) X 10"^°. (78) 

We have shown here SM predictions for these observables that do not include the 
correction r{ys). As r{yd) = 1 to an excellent approximation, the result for — )■ fi'^fi' 
can be directly compared with experiment. In order to obtain these results we have 
used the parametric formulae of 42 and updated the lattice QCD values of Fb^ ^ 43 
and the life-times tb 



44 



17 



. Details can be found in 
If the correction factor r(?/s) is taken into account the SM result in (77) changes to [17] 

(3.56 ±0.18) ■ 10~^ 



l3{Bs ^ fi+fi- 



)SM 



(79) 



It is this branching that should be compared in such a case with the results of LHCb 
given above. For the latest discussions of these issues see |l0|[ll|[37||42] . As discussed 
in 42 ,45 complete NLO electroweak corrections are still missing in this estimate. This 
result should be available in the near futureE] 



In our numerical results we will use B{Bs — )■ ) in (74) with B{Bs — ?■ ) given 



by (75) and r{ys) by (73) with also affected by NP effects. 
Combining the experimental and theoretical results quoted above gives 



R 



LHCb 



0.90l°j^ e [0.30, 1.80] (95% C.L). 



This range should be compared with its SM value, corresponding to P = 1, S" 
0r = 0: 

-RsM = 1. 



(80) 
and 

(81) 



4.5 Scenarios for S and P 

Clearly, the outcome of the results for the observables in question depends on the values 
of the muon couplings and whether a scalar or pseudoscalar boson is involved. Moreover, 
as we stressed already in Section [2] the exchanged mass eigenstate does not have to be a 



CP eigenstate and can have both scalar and pseudoscalar couplings to leptons. In 17 
a detailed classification of various possibilities beyond the dynamical model considered 
here has been made and the related purely phenomenological numerical analysis has 
been performed. Here we will make a classification that is particularly suited for the 
dynamical model considered by us. 



^Martin Gorbahn, private communication. 
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Pseudoscalar Scenario: 

In this scenario S = and P can be arbitrary complex number. We find then 



The branching ratio observable is given by 

'l + Vs cos(2v?p - 2ipBs 



R 



12) 



^3) 



This scenario corresponds to scenario A in p!?]. 
Scalar Scenario: 

In this scenario P = 1 and S can be arbitrary complex number. We find then 

COs2v9b^ - |S'pCOs(2v95 - 2ipBs) 



AT 



s 



lip. 



R 



1 + \S\^ 

-sin2(pBs - \S\'^sm{2(ps -2!fB, 
1 + \S\^ 

l + ys cos 2ipB, 



1 + 2/s 



+ \S\^ 



l-Vs cos(2ips - 2ipB, 



This scenario corresponds to scenario B in 17 
Mixed Scenario: 

We will consider a scenarios in which P is modified from its SM value and S is non-zero. 
As we want to discuss the case of a single new particle with spin 0, this means that this 
particle has both scalar and pseudoscalar couplings to muons. 

A simple scenario with both a scalar {H^) and pseudoscalar {A^) with approximately 
the same mass that couple equally to quarks and leptons up to the usual i factor in the 



pseudoscalar coupling has been recently considered in 17 . This scenario can be realized 



as a special limit in models like 2HDM and the MSSM with interesting consequences 
for ^^p. We refer to |17] for details. 
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5.1 Effective Hamiltonian for d 



For the study of Kl — )■ /i^yU and Kl — >■ Ti^ii decays we will need the relevant effective 



Hamiltonian. It can be obtained from the formulae of subsection 4.2 For completeness 
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we list here explicit formulae for operators and Wilson coefficients: 

Qg = {s^f.PLd){iY(), Q9 = {si^PRd){iri), (85a) 

Qio = {s-f^PLd){iYl5i), Q'lo = {s-f,PRd){iYl5i), (85b) 

Qs = m,{sPLd){U), Q's = m,{sPRd){ii), (85c) 

Qp = m,{sPLd){i^,i), Q'p = m,{sPRd){hd). (85d) 

Note that because of the sd ordering instead of qh scalar operators have L and R 
interchanged with respect to 6 — )■ s, (i transitions. 

The Wilson coefficients Cg and Cio do not receive any new contributions from scalar 



exchange and take SM values as given in (55). However, in order to include charm 



component in Kl we make replacement: 



ilvY^ixt) VvYoixt) + B^>NNL (86) 



where at NNLO 46 



>NNL = X'Pc{Y), Pc{Y) = 0.113 ± 0.017 . (87) 
The coefficients of scalar operators are: 

rris sin^ OwCs 

rus sm Uw(~'s - ^J72 TT^TT ' (89) 

sin^ 9^Cp=-^-^ 'v*V ' (^0) 
2^ ^, 11 A'^(H)Af(H) , , 

5.2 fi+fi- 

Only the so-called short distance (SD) part to a dispersive contribution to — )• 
can be reliably calculated. Therefore in what follows this decay will be treated only 
as an additional constraint to be sure that the rough upper bound given below is not 
violated. 



1 


1 


At{H)AfiH) 








1 


1 


A^i{H)Af{H) 


9sM 


Ml 


V^sVtd 


1 


1 


Af{H)Af{H) 


9sM 


Ml 




1 


1 


Ali{H)Af{H) 


9sM 


Ml 


Vt*sVtd 
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The relevant branching ratio can be obtained by first introducing: 



P{K) = Cio - C[ 



10 



2m„ nid + 



iCp-C'p) 



SiK) 



mj^ 2m^ ma + rris 



We then find 



B{Kl )sD 



47^5 



4m?, 



^ sin^ d 



m 



w 



(92) 
(93) 

(94) 



K 



X 



3« V*VuP 



+ 



and 3? ^ 3 for Ks — t- fi^fi decay. We recall that Ciq does not receives any contribution 



from scalar exchanges and includes also SM charm contribution as given in (86). C[q = 
for scalar exchanges. 

Equivalently we can write 



B{Kl )sD 



where 



a^B{K+ fi+u) t{Kl) 



Q f V*VuS 



(95) 



(96) 



A%2 r{K+)' 

The extraction of the short distance part from the data is subject to considerable 



uncertainties. The most recent estimate gives 47 



B{Kl ^ /x+yU-)sD < 2.5 ■ 10 



-9 



(97) 



to be compared with (0.8 ± 0.1) ■ 10"^ in the SM 46 . 



5.3 Kl 7r^£+i- 

The rare decays Ki — )■ 7r°e~*"e~ and Ki — )■ 7r°yU+/i^ are dominated by CP-violating 
contributions. The indirect CP-violating contributions are determined by the measured 
decays Ks — J- 7r°£+£~ and the parameter ek in a model independent manner. It is the 
dominant contribution within the SM where one finds |48] 

B{Kl ^ 7r°e+e-)sM = 3.54;°|^ (1.56;°;^^) . ^q-ii ^ (gg) 
BiK, ^ 7rW-)sM = 1.41^°:i (0.95^°|?) • lO-^^ , (99) 
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with the values in parentheses corresponding to the destructive interference between 
directly and indirectly CP- violating contributions. The last discussion of the theoretical 
status of this interference sign can be found in 49 where the results of 50-52 are 
critically analysed. From this discussion, constructive interference seems to be favoured 
though more work is necessary. In spite of significant uncertainties in the SM prediction 
we will investigate how large the scalar contributions to these decays are still allowed 
by present constraints. To this end we will confine our analysis to the case of the 
constructive interference between the directly and indirectly CP-violating contributions. 

The present experimental bounds 



-B(irL^7r°e+e")exp<28-10"^i 53, -6(7^^ ^ 7rVV~)cxp < 38 ■ 10"^^ |54 



(100) 

are still by one order of magnitude larger than the SM predictions, leaving thereby 
large room for NP contributions. While in the case of Z' models large enhancements 
of branching ratios were not possible due to constraints from data on — )■ vr+z/i/ [I], 
this constraint is absent in the case of scalar contributions and it is of interest to see 
by how much the branching ratios can be enhanced in the models considered here still 



being consistent with all data, in particular with the bound in (97) 



In the LHT model the branching ratios for both decays can be enhanced at most by a 
factor of 1.5 ^[S]. Slightly larger effects are still allowed in Randall-Sundrum models 
with custodial protection (RSc) for left-handed couplings (6]. Even larger effects are 



found if the custodial protection is absent 55 



Probably the most extensive model independent analysis of decays in question has 



been performed in 48 , where formulae for branching ratios for both decays in the 
presence of new operators have been presented. These formulae have been already 
used in [4 56 for the LHT model and in [6] in the case of RSc. In the LHT model. 



where only SM operators are present the effects of NP can be compactly summarized 
by generalization of the real SM functions lo(xt) and Zq^Xi) to two complex functions 
Yk and Zk-, respectively. As demonstrated in the context of the corresponding analysis 
within RSc |6] , also in the presence of RH currents two complex functions Yk and Zk 
are sufficient to describe jointly the SM and NP contributions. Consequently the LHT 
formulae (8.1)-(8.8) of |4| with Yk and Zk given in (88) and (89) of (l) can be used 
in the context of tree-level gauge boson exchanges. The original papers behind these 
formulae can be found in |48l|50l[5Tll57l|58l . 



The case of scalar contributions is more involved. In order to use the formulae of 48 
for scalar contributions we introduce the following quantities: 

_ 1 r^yY,{x,) 
'"'^--2^ sin^e^ 1.4-10-^' 

yp + y'p - ys + y's /mo^ 

yp = t; , ys = — ^ — (102) 
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with Hi related to the Wilson coefficients in the present paper as follows: 

M^sin^ M^sm^ ^^^^^ 

mi mi 

with analogous formulae for primed coefficients. Here mi stands for mg and as the 



authors of 48 anticipating hehcity suppression included these masses already in the 



effective Hamiltonian. 



Using 48 we find then corrections from tree-level and exchanges to the branching 



ratios that should be added directly to SM results in (98) and (99): 



AB'p''- = (1.9 u-TA^iyp) + 0.038 myp))') ■ lO''', (104) 

AB^p"' = (0.26 ujA'^iyp) + 0.0085 {'^{yp))^) ■ 10"^^ (105) 

ABf^' = (1.5 ^(ys) + 0.0039 mys))') ■ W'^^ (106) 

ABf^~ = (0.04 ^{ys) + 0.0041 {?fi{ys)Y) ■ 10"''. (107) 
Note that in the absence of helicity suppression the large suppression factors above are 



canceled by the conversion factors in (103). 



The numerical results for these new contributions are given in Section [9] 

6 General Structure of New Physics Contributions 
6.1 Preliminaries 

We have seen in Section |2] that the small number of free parameters in each of LHS, 
RHS, LRS and ALRS scenarios allows to expect definite correlations between flavour 
observables in each step of the strategy outlined there. These expectations will be 
con&med through the numerical analysis below but it is instructive to develop ffist a 
qualitative general view on NP contributions in different scenarios before entering the 
details. 

First, it should be realized that the confrontation of correlations in question with future 
precise data will not only depend on the size of theoretical, parametric and experimental 
uncertainties, but also in an important manner on the size of allowed deviations from SM 
expectations. The latter deviations are presently constrained dominantly by AF = 2 



observables and B — )■ Xs'j decay. But as already demonstrated in 1 7 8,12 with the 
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the new data from the LHCb, ATLAS and CMS at hand also the decays Bg^d fJ'~^fJ'~ 
and b — )■ si^i~ begin to play important roles in this context. We will see their impact 
on our analysis as well. 

Now, in general NP scenarios in which there are many free parameters, it is possible with 
the help of some amount of fine-tuning to satisfy constraints from AF = 2 processes 
without a large impact on the size of NP contributions to AF = 1 processes. However, 
in the case at hand in which NP in both AF = 2 and AF = 1 processes is governed by 
flavour changing tree-diagrams, the situation is different. Indeed, due to the property 
of factorization of decay amplitudes into vertices and the propagator at the tree-level, 
the same quark flavour violating couplings and the same mass enter AF = 2 and 
AF = 1 processes undisturbed by the presence of fermions entering the usual box and 
penguin diagrams. Let us exhibit these correlations in explicit terms. 



6.2 AF = 1 vs. AF = 2 Correlations 

In order to obtain transparent expressions we rewrite various contributions [AS{K)]xb 
and [AS'(-Bg)]AB with A, B = L, R to AF = 2 amplitudes as follows 

[A5W]ab = -j^^ (108) 

[A5(i?,)]AB - (109) 



where the quantities r (M) can be found by comparing these expressions with (33) 



(34), (40) and (41) and analogous expressions for the contributions of operators Q'l 
They depend on low energy parameters, in particular on the meson system and loga- 
rithmically on Mh- The latter dependence can be neglected for all practical purposes 
as long as Mh is above several hundreds of GeV and still in the reach of the LHC. We 
collect the values of r^^(M) in Table [3] 

Defining 



Cfp = rui sin^ OwC^^p (1 10) 



we can then derive the following relations between the Wilson coefficients Cg^p entering 
the AF = 1 processes and the shifts [AS'(M)]ab in AF = 2 processes which are inde- 
pendent of any parameters like but depend sensitively on Mh and on the couplings 
Ag^p^H^j^ In particular they do not depend explicitly on whether SI or S2 scenarios 
for |V^fe| defined below are considered. This dependence is hidden in the allowed shifts 
in [AS{K)]ab and [AS{Bd)]AB both in magnitudes and phases. We have theiij^ 

^Similar relations have been derived in jlj in the context of Z' models. 

^The numerical values on the r.h.s of these equations correspond to Mh — 1 TeV. 
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LL/RR 


LR 


K 


960 


-5700 


Bd 


51 


-310 


Bs 


50 


-300 



Table 3: r^^(M) in units of TeV^ as defined m Eqs. (w^ and |7og| j for Mh = 1 TeV. 



= 0.18A^|^p(i/), 
-0.005[A^J,(i/)]^ 



For Bn we have to make the following replacements in the formulae above: 



[AS{K)Ub [^SiB,)U, 



Cs,p{Bq 



^[AS{B,)]l^ MHgl^./r^AB7) 
CsAB,)C'sAB,) [A^sAH)f 



[^S{B,)]l^ 



MWsur'^HB,) 



: Q.l^AfAH), 
O-YSA^sAH), 
-0.1[A^^J,(/7)]^ 



6.3 Implications 



(111) 
(112) 
(113) 

(114) 
(115) 

(116) 

(117) 



Inspecting these formulae we observe that if the SM prediction for Ek is very close to its 
experimental value, AS{K) cannot be large and consequently at first sight the values of 
the Wilson coefficients Agp{K) cannot be large implying suppressed NP contributions 
to rare K decays unless H couplings to charged leptons in the final state are enhanced, 
although this enhancement can be bounded by rare Bg^d decays. Further details depend 
on the value of Mh- While in Z' models the present theoretical and parametric uncer- 
tainties in ek and AMk still allow for large effects in rare K decays both in SI and S2 
scenarios, this turns out not to be the case in the models considered here. 



Similarly in the Bd and Bg systems if the SM predictions for AMg^d^ S^Ks ^'^^ 



■>ip(f, 



are very close to the data, it is unlikely that large NP contributions to rare Bd and Bs 
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decays, in particular the asymmetries >S'^+^_ , will be found, unless again H couplings to 
charged leptons in the final state are enhanced. Here the situation concerning theoretical 
and parametric uncertainties is better than in the K system and the presence of several 
additional constraints from 6 — )■ s transitions allows to reach in the Bg system clear cut 
conclusions. 

In this context it is fortunate that within the SM there appears to be a tension between 
the values of Ek and S^Ks 33 , 59 so that some action from NP is required. Moreover, 



parallel to this tension, the values of \Vub\ extracted from inclusive and exclusive decays 



differ significantly from each other. For a recent review see 60 . 

If one does not average the inclusive and exclusive values of \Vub\ and takes into account 
the tensions mentioned above, one is lead naturally to two scenarios for NP: 

• Exclusive (small) \Vub\ Scenario 1: {ekI is smaller than its experimental de- 
termination, while S^Kg is rather close to the central experimental value. 

• Inclusive (large) \Vub\ Scenario 2: \ek\ is consistent with its experimental 
determination, while S^Ks is significantly higher than its experimental value. 

Thus depending on which scenario is considered, we need either constructive NP con- 
tributions to l^i^l (Scenario 1) or destructive NP contributions to S^pxs (Scenario 2). 
However this NP should not spoil the agreement with the data for S^Ks (Scenario 1) 
and for \ek\ (Scenario 2). 

While introducing these two scenarios, we should emphasize the following difference 
between them. In Scenario 1, the central value of \ek\ is visibly smaller than the very 
precise data but the still significant parametric uncertainty due to \Vcb\^ dependence in 
l^j^l and a large uncertainty in the charm contribution found at the NNLO level in 32 
does not make this problem as pronounced as this is the case of Scenario 2, where large 
\Vub\ implies definitely a value of S^Ks that is by 3cr above the data. 

Our previous discussion allows to expect larger NP effects in rare Bd decays in scenario 
S2 than in SI. This will be indeed confirmed by our numerical analysis. In the K 
system one would expect larger NP effects in scenario SI than S2 but the present 
uncertainties in Ek and AM^ do not allow to see this clearly. The Bg system is not 
affected by the choice of these scenarios and in fact our results in SI and S2 are basically 
indistinguishable from each other as long as there is no correlation with the Bd system. 
However, we will demonstrate that the imposition of t/(2)^ symmetry on H couplings 
will introduce such correlation with interesting implications for the Bg system. 

We do not include B ^ r^Ur in this discussion as NP related to this decay has nothing 
to do with neutral scalars, at least at the tree-level. Moreover, the disagreement of 
the data with the SM in this case softened significantly with the new data from Belle 
Collaboration IgTI. The new world average provided by the UTfit collaboration of 



B{B+ r+z/)exp = (0.99 ± 0.25) x lO""^ 62 is in perfect agreement with the SM in 



scenario S2 and only by 1.5cr above the SM value in scenario SI. 
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Evidently \Vub\ could be some average between the inclusive and exclusive values, in 
which significant NP effects will be in principle allowed simultaneously in K and 
decays. This is in fact necessary in NP scenarios in which NP effects to AF = 2 
processes are negligible and some optimal value for \Vub\, like 0.0037 is chosen in order 
to obtain rough agreement with the data. But then one should hope that future data, 
while selecting this value of \Vub\, will also appropriately imply a higher experimental 
value of S^Ks ci-nd new lattice results will bring modified non-perturbative parameters in 
the remaining AF = 2 observables so that everything works. This is the case of a recent 



analysis of FCNC processes within a model for quark masses 63 . This discussion shows 



importantance of the determination of the value of \Vub\ and of the non-perturbative 
parameters in question (see article by A. Buras in |64J). 

As already remarked above, the case of Bs mesons is different as the B^ — B^ system 
is not involved in the tensions discussed above. Here the visible deviation of the AM^ 
in the SM from the data and the asymmetry S^^, still being not accurately measured, 
govern the possible size of NP contributions in rare decays. 



6.4 Dependence on Mh 



The correlations between AF = 1 and AF = 2 derived in subsection |6.2| imply that 
when free NP parameters have been bounded by AF = 2 constraints, the Wilson 
coefficients of scalar operators are inversely proportional to M^- This means that 
in the case of NP contributions significantly smaller than the SM contributions in P, 
the modifications of rare decay branching ratios due to NP will be governed by the 
interference of SM and NP contributions. Consequently such contributions to branching 
ratios will also be inversely proportional to M^- If NP contribution to P is of the size of 
the SM contributions than this law will be modified and NP contributions will decrease 
faster with increasing Mh- On the other hand in the absence of interference between 
NP and SM contributions, as is the case of S, the NP modifications of branching ratios 
will decrease as l/Mfj. Consequently, we expect that for sufficiently large Mh only 
NP contributions in P, as in Z' scenarios, will matter unless the scalar couplings are 
very much enhanced over pseudoscalar ones. Evidently, for low values of Mh the 5* 
contributions could be relevant. Here in principle a SM Higgs, being a scalar, could 
play a prominent role, but as we will demonstrate below this can only be the case for 
AF = 2 transitions. 

While Mh could still be as low as few hundreds of GeV, in order to cover a large set 
of models, we will choose as our nominal value Mh = 1 TeV. With the help of the 



formulae in subsection 6.2 it should be possible to estimate approximately, how our 
results would change for other values of Mh- In this context it should be noted that 
any change of Mh can be compensated by the change in couplings A^^p unless these 
couplings are predicted in a given model or are known from other measurements. 



With this general picture in mind we can now proceed to numerical analysis. 
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7 Strategy for Numerical Analysis 
7.1 Preliminaries 

Similarly to our analyses in |l(,|12j it is not the goal of the next section to present a full- 
fledged numerical analysis of all correlations including present theoretical, parametric 
and experimental uncertainties as this would only wash out the effects we want to 



emphasize. Yet, these uncertainties will be significantly reduced in the coming years 65 
66 and it is of interest to ask how the H scenarios considered here would face precision 
flavour data and the reduction of hadronic and CKM uncertainties. In this respect, as 
emphasized above, correlations between various observables are very important and we 
would like to exhibit these correlations by assuming reduced uncertainties in question. 

Therefore, in our numerical analysis we will choose as nominal values for three out of 
four CKM parameters: 



IK.I = 0.2252, 



IK 



cb\ 



0.0406, 



7 = 68°, 



(118) 



and instead of taking into account their uncertainties directly, we will take them effec- 
tively at a reduced level by increasing the experimental uncertainties in AMg^^ and Ek- 
Here the values for \Vus\ and \Vcb\ have been measured in tree level decays. The value 
for 7 is consistent with CKM fits and as the ratio AM^/AM^ in the SM agrees well 
with the data, this choice is a legitimate one. Other inputs are collected in Table |4j 
For \Vuh\ we will use as two values 



IK 



ub I 



3.1 ■ 10 



-3 



IK 



uh I 



4.0 ■ 10 



-3 



;ii9) 



that are in the ballpark of exclusive and inclusive determinations of this CKM element 
and representing thereby SI and S2 scenarios, respectively. 



Having fixed the three parameters of the CKM matrix to the values in (118), for a given 
\Vub\ the "true" values of the angle /3 and of the element |Kd| are obtained from the 
unitarity of the CKM matrix: 



IK 



td\ 



where 



iKisI |Kfe|-Rt5 



l + Rl- 2Rb cos 7 , 



cot [3 



— Rb cos 7 
Rb sin 7 



Rb 



2 



1 


|K,| 


A 


\Vcb\ 



(120) 



(121) 



In Table |5| we summarize for completeness the SM results for \eK 
and B{B 
for 



AM. 



s4, ^sin2/3),^^^ 
and choosing the two values 



obtained from (120), setting 7 = 68 
iKifel in ( |119 ). We observe that for both choices of \Vub\ the data show significant 
deviations from the SM predictions but the character of the NP which could cure these 
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Gf = 1.16637(1) X 10- 
Mw = 80.385(15) GeV 
sin^Ow = 0.23116(13) 
a{Mz) = 1/127.9 
a,(Mz) = 0.1184(7) 



GeV 



=2" 



m„(2GeV) = (2.1 ±0.1)MeV 
md{2 GeV) = (4.73 ± 0.12) MeV 

m,(2 GeV) = (93.4 ± 1.1) MeV 



(1.279 ±0.013) GeV 

+0.18 , 



4. 19--- GeV 



mb{mb) = 
mt{mt) = 163(1) GeV 
Mt = 173.2 ± 0.9 GeV 
rriK = 497.614(24) MeV 
Fk = 156.1(11) MeV 
Bk = 0.767(10) 

= 0.94(2) 
r/i = 1.87(76) 

= 0.5765(65) 
r^s = 0.496(47) 

AMk = 0.5292(9) x lO^^pg-i 
IekI = 2.228(11) X IQ-^ 



68 



33 



B{B 



(3.55 ±0.24 ±0.09) x 
= (0.99 ±0.25) X 10"^ 



10" 



tb± = (1641 ± 8) X lO^ps 



67 
67 
67 

E 

67 



68 
68 

68 



69 

S 

70 

n 
m 

68 
68 
34 
32 
29 

H 
67 



67 
62 



44 



ruB, 

Bb. 
Bb. 

Fb. 



5279.5(3) MeV 
= 5366.3(6) MeV 
= (188 ±4) MeV 
= (225 ± 3) MeV 
= 1.26(11) 
= 1.33(6) 
/Bb^= 1.05(7) 

Bb^ = 226(13) MeV 

Fbs\/Bbs = 279(13) MeV 

^ = 1.237(32) 

riB = 0.55(1) 

AMd = 0.507(4) ps"^ 

AM, = 17.72(4) ps-i 

S^Ks = 0.679(20) 

= 0.001 ± 0.100 
Ar, = 0.116 ±0.019 
t{Bs) = 1.503(10) ps 
r(5,) = 1.519(7) ps 



29 



67 

67' 
43' 

68 
68 
68 



68 



68 
W 
30 

34 
44 
Wi 
24 



24 



\Vus\ = 

\Vcb\ = 
I "j/incl. I 
I ^ub I 

I T/OXCl. 



0.2252(9) |67| 
(40.9 ± 1.1) X 10-3 |72| 
= (4.41 ±0.31) X 10-3 p] 
= (3.23 ±0.31) X 10-fl 



Table 4: Values of the experimental and theoretical quantities used as input parameters. 





Scenario 1: 


Scenario 2: 


Experiment 




1.72(22) ■ 10-3 


2.15(32) ■ 10-3 


2.228(11) X 10-3 


(sin 2/3)tme 


0.623(25) 


0.770(23) 


0.679(20) 


AMs [ps-i] 


19.0(21) 


19.0(21) 


17.73(5) 


AMd [ps-i] 


0.56(6) 


0.56(6) 


0.507(4) 




0.62(14) ■ 10-^ 


1.02(20) ■ 10-^ 


0.99(25) X 10-4 



Table 5: SM prediction for various observables for \Vub\ 
and 7 = 68° compared to experiment. 



3.1-10-3 and IK 



ub I 



4.0-10- 



tensions depends on the choice of \Vub\ as already discussed in detail in 73 and in the 
previous section. 



What is striking in this table is that the predicted central values of AM, and AMd, 
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although slightly above the data, are both in good agreement with the latter when 
hadronic uncertainties are taken into account. In particular the central value of the 
ratio AMs/AMrf is very close to the data: 

-— i =34.5 ±3.0 exp: 35.0 ± 0.3 . 122 

These results depend on the lattice input and in the case of AM^ on the value of 7. 
Therefore to get a better insight both lattice input and the tree level determination of 
7 have to improve. 

Similarly to the anatomy of Z' models in we will deal with two scenarios for and 
four scenarios LHS, RHS, LRS, ALRS for flavour violating couplings of H to quarks. 
Thus for a given scalar or pseudoscalar we will deal with eight scenarios of flavour 
violating if-physics to be denoted by 

LHSl, LHS2, RHSl, RHS2, LRSl, LRS2, ALRSl, ALRS2 (123) 

with SI and S2 indicating the \yuh\ scenarios. With the help of scalar, pseudoscalar 



and mixed scenarios for leptonic couplings introduced in Subsection |4.5| in each case, 
we will be able to get the full picture of various possibilities. 



We should emphasize that in each of the scenarios listed in (123), except for leptonic 
couplings, we have only two free parameters describing the if-quark couplings in each 
meson system except for the universal Mh- Therefore, as in the case of Z' models it is 
possible to determine these couplings from flavour observables (see Section [2]) provided 
flavour conserving H couplings to muons and Mu are known. While in the SM and some 
specific models scalar couplings are known, in the present analysis we want to be more 
model independent. While we will get some insight about them from — ?■ yu^/i", 
determining them in purely leptonic processes increases the predictive power of the 
theory. 

Following Step 2 of our general strategy of Section [2} in what follows we will assume 
that Ap^(if) and Ag'^lH) have been determined in purely leptonic processes. For 
definiteness we set the lepton couplings at the following values 

Af{H) = ±0.020, Af{H) = 0.040. (124) 

As we will demonstrate in the course of our presentation these values are consistent with 
the allowed range for B{Bs — ?■ yU'^'Ai") when the constraints on the quark couplings from 
— B^ are taken into account and Mh = 1 TeV. The reason for choosing the scalar 
couplings to be larger than the pseudoscalar ones is that they are weaker constrained 
than the latter because the scalar contributions do not interfere with SM contributions. 
Note that because of the lack of this interference, the values of S are simply proportional 
to A'g^{H) and it is straightforward to obtain S contributions for different values of this 
coupling. 
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These couplings should be compared with SM Higgs couplings 

[Af (if)]sM = 0, [Af (i/)]sM = 1.2 X 10-3. (125) 

The smallness of these couplings precludes any visible SM Higgs effects in B{Bs — t- 
even in the presence of flavour violating quark couplings. Indeed the latter have 
to be smaller than in the case of heavy scalars, suppressing NP contributions to the 
branching ratios even further. 



Concerning the signs in (124), the one of Ag^{H) is irrelevant as only the square of this 
coupling enters various observables. The sign of A^p{H) has an impact on the interfer- 
ence of pseudoscalar and SM contributions and is thereby crucial for the identification 
of various enhancements and suppressions with respect to SM branching ratios and CP 
asymmetries. Consequently it plays a role of our search for successful oases in the space 
of parameters. 



7.2 Simplified Analysis 

As in [I] we will perform a simplified analysis of Ek, AM^^^, S^Ks and S^^ in order to 
identify oases in the space of new parameters (see Section [2]) for which these five ob- 
servables are consistent with experiment. To this end we set all other input parameters 
at their central values but in order to take partially hadronic and experimental uncer- 



tainties into account we require the theory in each of the eight scenarios in (123) to 
reproduce the data for sk within ±10%, AM^^^ within ±5% and the data on Sr^pKs and 
S^^ within experimental 2a. We choose larger uncertainty for ek than AM^^rf because 
of its strong \Vcb\^ dependence. For AMk we will only require the agreement within 
±25% because of potential long distance uncertainties. 

Specifically, our search is governed by the following allowed range^ 

16.9/ps < AM, < 18.7/ps, -0.20 < S^^ < 0.20, (126) 
0.48/ps < AMd < 0.53/ps, 0.64 < S^Ks < 0-72. (127) 

0.75 < ^ < 1.25, 2.0 X 10-3 < \ek\ < 2.5 x lO-^. (128) 



The search for these oases in each of the scenarios in (123) is simplified by the fact 
that for fixed Mh each of the pairs (AM,, S^^), (AM^^, S^Ks) and (AM^^ , \ek\) depend 
only on two variables. The fact that in the K system we have only one powerful 
constraint at present is rather unfortunate. Moreover, in the models considered the 
decays — > -k^vv and Kl — )■ ti^vv cannot help unless charged Higgs contributions 



■^When using the constraint from we take into account that only mixing phase close to its 
SM value is allowed thereby removing some discrete ambiguities. The same is done for S^pKs- 
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are considered, which is beyond the scope of the present paper. While the constraint 



(97) on Kl — )■ could have in principle an impact on our search for oases, we have 



checked that this is not the case. 

In what follows we will first for each scenario identify the allowed oases. As in the case of 



Z' models there will be in principle four oases allowed by the constraints in (126)- (127). 
However, when one takes into account that the data imply the phases in S^^ and S^Ks 
to be close to the SM phases, only two big oases are left in each case. Similarly the sign 
of Ek selects two allowed oases. In order to identify the final oasis we will have to invoke 
other observables, which are experimentally only weakly bounded at present. Yet, our 
plots will show that once these observables will be measured precisely one day not only 
a unique oasis in the parameter space will be identified but the specific correlations in 
this oasis will provide a powerful test of the flavour violating H scenarios. 



8 An Excursion through H Scenarios 

8.1 The LHSl and LHS2 Scenarios 
8.1.1 The Bs Meson System 

We begin the search for the oases with the Bg system as here the choice of \Vub\ is 
immaterial and the results for LHSl and LHS2 scenarios are almost identical. Basically 
only the asymmetry S^^ within the SM and \Vts\ are slightly modified because of the 
unitarity of the CKM matrix. But this changes S^^ in the SM from 0.032 to 0.042 and 
can be neglected. 

The result of this search for Mh = 1 TeV is shown in Fig. [2| where we show the allowed 
ranges for (523,523). The red regions correspond to the allowed ranges for AMg, while 
the blue ones to the corresponding ranges for S^^p. The overlap between red and blue 
regions identifies the oases we were looking for. We observe that the requirement of 
suppression of AMg implies S23 7^ 0. 

Comparing Fig. [2] with the corresponding Z' result of Fig. 2 in fl] we observe that 
the phase structure is identical to the one found in the case of Z' but the values of 
S23 are smaller. This behaviour is easy to understand. While the tree diagram with 
scalar exchange has the overall sign opposite to the one of a gauge boson exchange, this 
difference is canceled by the opposite signs of the matrix element of the leading operator 
Qf^^ and Q^^^ in the case of H and Z' exchange, respectively. But the absolute value 
of (Qi^^) is larger than of (QY^^) and consequently S23 in the Higgs case has to be 
smaller than in the Z' case in order to fit data. We find that this suppression of S23 
that enters quadratically in AM,, amounts roughly to a factor of 1.5. 

In view of this simple change we do not show the table for the allowed ranges for ^23 
and §23- They are obtained from Table 5 in ^ by leaving ^23 unchanged and rescaling 
§23 down by a factor of 1.5. 
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AM, S S,.,,>, LHS1 




0.000 0.002 0.004 0.006 O.OOS 0.010 0.012 0.014 



Figure 2: Ranges for AMg (red region) and S^rp (blue region) for Mh = 1 TeV in LHSl 



satisfying the hounds in Eq. (126). 



Inspecting Fig. |2] we observe the following pattern: 

• For each oasis with a given ^23 there is another oasis with ^23 shifted by 180° but 
the range for §23 is unchanged. This discrete ambiguity results from the fact that 
AMs and S^^ are governed by 2523. However, as we will see below this ambiguity 
can be resolved by other observables. Without the additional information on 



phases mentioned in connection with constraints (126)-(128) one would find two 
additional small oases corresponding roughly to NP contribution to twice as 
large as the SM one but carrying opposite sign. But taking these constraints on 
the phases into account removes these oases from our analysis. 

• The increase of Mh by a given factor allows to increase §23 by the same factor. 
This structure is evident from the formulae for AS{Bs). However, the inspection 
of the formulae for AF = 1 transitions shows that this change will have impact 
on rare decays, making the NP effects in them with increased Mh smaller. This 
is evident from the correlations derived in Subsection 16.21 

We will next confine our numerical analysis to these oases, investigating whether some 
of them can be excluded by other constraints and studying correlations between various 
observables. To this end we consider in parallel pseudoscalar and scalar scenarios setting 



the lepton couplings as given in (124). In addition to the general case corresponding 
to the oases just discussed we will present in plots the results obtained when the f/(2)^ 
symmetry is imposed on Bg and B^i systems. This case will be discussed in detail at 
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the end of this subsection but to avoid too many plots and to show the impact of 
this symmetry we will already include the results in discussing the results without this 
symmetry. Our colour coding will be as follows: 



• In the general case blue and purple allowed regions correspond to oases with small 
and large 623, respectively. However, one should keep in mind the next comment. 

• In the U(2)^ symmetry case, the allowed region will be in magenta and and cyan 
for LHSl and LIIS2, respectively, as in this case even in the Bs system there is 
dependence on \Vub\ scenario. These regions are subregions of the general blue or 
purple regions so that they cover some parts of them. 



In order to justify the values for the leptonic couplings in 
function of A'^ and A'^^ in LHSl for the pseudoscalar and scalar scenario, respectively. 
In ALRS effects are smaller and in LRS R does not depend on Ap^ and Ag'^. We 
observe that for equal scalar and pseudoscalar couplings, the effects are significantly 
larger in the A'^ case and this is the reason why we have chosen the scalar couplings to 
be larger. 

There are two striking differences between and cases originating in the fact that 
pseudoscalar contributions interfere with the SM contribution, while this is not the case 
for a scalar: 



124) we show in Fig. \3\R as 



• While in the H'^ case R can only be enhanced, it can also be suppressed in the 
A^ case. This difference could play an important role one day. 

• In the A° case the result depends on the oasis considered and the sign of Ap. 
However changing simultaneously the sign of Ap and the oasis leaves R invariant. 
In the case R is independent of the oasis considered and of the sign of A5. 

In Fig. |4] (left) we show vs S^p^ in the A^ case. In the same figure (right) we 

show the correlation between B{Bs — )■ /i"*"/^") and -S'^^j^ We observe that for largest 
allowed values of S^^ the asymmetry can be as large as ±0.5. Also the effects in 

B{Bs — 7- /i^yu^) are expected to be sizable for the chosen value of muon coupling. 

Comparing the plots in Fig. |4] with the corresponding results for Z' in Fig. 3 of [l] we 
observe striking differences which allow to distinguish the case of tree-level pseudoscalar 
exchange from the heavy gauge boson exchange: 

• In the A^ case the asymmetry >S'^+^_ can be zero while this was not the case in 
the Z' case where the requirement of suppression of AMg directly translated in 

^The central values for B{Bd M+//~)sm = 1-0 x 10~^° and B{Bs H+f^^)sM = 3.45 x IQ-^ 
shown in the plots correspond to fixed CKM parameters chosen by us and differ from the ones listed 



in ( 78 1 and ( 79 ) but are fully consistent with them. 
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ICC 




Figure 3: R as a function of Ap^ and Ag'^ in LHSl. Left: A'^ case; right: case. 
Gray region: 95% CL of R. 



^li+ii- being non-zero. Consequently in the Z' case the sign of S^+^^ could be 
used to identify the right oasis. The left plot in Fig. |4] clearly shows that this is 
not possible in the case. We also find that while in the Z' case the asymmetry 
could reach values as high as ±0.9, in the A^ case |S'^+^_| can hardly be 
larger than 0.5. 

• On the other hand we observe that in the A^ case the measurement of B{Bs — )■ 
uniquely chooses the right oasis. The enhancement of this branching ratio 
relatively to the SM chooses the blue oasis while the suppression the purple one. 
This was not possible in the Z' case. The maximal enhancements and suppressions 
are comparable in both cases but finding B{Bs — )■ fi'^fi~) close to SM value would 
require in the A^ case either larger or smaller muon coupling. 

We observe that the roles of and B{Bs — > fi^fJ^') in searching for optimal oasis 

have been interchanged when going from the Z' case to the A° case. While B{Bs 
IJL^ identifies the oasis the correlation of S^+^. vs. S^^ constitutes an important test 
of the model. While in the blue oasis S^^j, increases (decreases) uniquely with increasing 
(decreasing) 5*^+^- , in the purple oasis, the increase of 5*^^ implies uniquely a decrease 
of 'S'^+^_. Therefore, while S^+^^ alone cannot uniquely determine the optimal oasis, it 
can do in collaboration with S^p^. 

If the favoured oasis will be found to differ from the one found by means of B{Bs — t- 
H'^ one day the model in question will be in trouble. Indeed, let us assume that 
B{Bs — )■ fi^fi^) will be found above its SM value selecting thereby blue oasis. Then the 
measurement of S^^p will uniquely predict the sign of Moreover, in the case of 

S^^ sufficiently different from zero, we will be able to determine not only the sign but 
also the magnitude of S'f,+ „_. 

These striking differences between the A^-scenario and Z'-scenario can be traced back 
to the difference between the phase of the NP correction to P in these two NP scenarios. 
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Figure 4: >S'^+^_ versus S^^ (left) and S^^ versus B{Bs — )■ yU^yU ) (right) for Mh 



1 TeV in LHS for two oases as explained in the text and A case. The blue and purple 
regions are almost identical for LHSl and LHS2. The magenta region corresponds to 
the f/(2)3 limit for LHSl and the cyan region for LHS2 (see Sect. \8.1.^ . Gray region: 
exp la range B{Bs 



(3.2 



+1.5\ 
1.2) 



10 . Red point: SM central value. 



As the oasis structure as far as the phase ^23 is concerned is the same in both scenarios 
the difference enters through the muon couphngs which are imaginary in the case of 
-scenario but real in the case of Z' . This is in fact the main reason why the structure 
of correlations in both scenarios is so different. Taking in addition into account the sign 
difference between Z' and pseudoscalar propagator in the the h — > sfi~^fi~ amphtude, 
which is now not compensated by a hadronic matrix element, we find that 

P{Z') = 1 + rz'c'^^', P{A^) = 1 + TAoe'^^'' (129) 

with 

TT 

rz' ^ tao, Sz' = S23 - 13s, Sao = 6z' - -■ (130) 

Therefore with ^23 of Fig. [2] the phase dz' is around 90° and 270° for the blue and purple 
oasis, respectively. Correspondingly 5^0 is around 0° and 180*^. This difference in the 
phases is at the origin of the differences listed above. In particular, we understand now 
why the CP asymmetry >S'*+^_ can vanish in the A^ case, while it was always different 
from zero in the Z'-case. What is interesting is that this difference is just related to the 
different particle exchanged: gauge boson and pseudoscalar. We summarize the ranges 
of 5z' and in Table |6| 

The power of the correlations in question in distinguishing between various scenarios 
is further demonstrated when we consider the case of a scalar in which there is no 
interference with the SM contribution. In Fig. [5] we show the corresponding results in 
the if" case. We observe the following differences with respect to Fig. |4j 

• B{Bs — )■ yU^/i^) can only be enhanced in this scenario and this result is inde- 
pendent of the oasis considered. Thus finding this branching ratio below its SM 
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Oasis 


6z' 


Sao 


Bs (blue) 
Bs (purple) 


50° - 130° 
230° - 310° 


-40° - (+40°) 
140° - 220° 


Bd (SI) (yellow) 
5, (SI) (green) 

Bd (S2) (yellow) 
Bd (S2) (green) 


57° - 86° 
237° - 266° 
103° - 125° 
283° - 305° 


-33° - (+4°) 
147° - 176° 

13° - 35° 
193° -215° 


U{2f (SI) (blue, magenta) 
f/(2)^ (SI) (purple, magenta) 
[7(2)3 (S2) (blue, cyan) 
U{2f (S2) (purple, cyan) 


55° - 84° 
235° - 264° 
101° - 121° 
291° -301° 


-35° - (-6°) 
145° - 174° 

11° -31° 
201° - 211° 



Table 6: Ranges for the values of 6z' and 5ao as defined in (129) for the Bg and Bd 
systems and various cases discussed in the text. Also the result for U{2)^ models is 
shown. 
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Figure 5: 5'^+^- versus S^^ (left) and S^^j, versus B{Bs jJi^jJi ) (right) for = 
1 TeV in LHSl for two oases as explained in the text and case (the two oases overlap 
here). The magenta region corresponds to the U{2)^ limit for LHSl and the cyan region 
for LHS2. Gray region: exp la range B{Bs — )• /i^/i^) 
central value. 



{?,.2t\l) ■ 10-^ Red point: SM 



value would favour the pseudoscalar scenario over scalar one. But the enhance- 
ment is not as pronounced as in the pseudoscalar case because the correction to 
the branching ratio is governed here by the square of the muon coupling while in 
the pseudoscalar case the correction was proportional to this coupling due to the 
interference with the SM contribution which is absent here. 

• Concerning CP-asymmetries similarly to the branching ratio there is no depen- 
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LHS & 4° LHS & H° 




Figure 6: A'^j- versus S^fj, for = 1 TeV in LHSl, blue and purple oases and 
case (left) and case (right). In case the two oases overlap. The magenta region 
corresponds to the U (2)^ limit for LHSl and the cyan region for LHS 2. Red point: SM 
central value. 



dence on the oasis considered but more importantly 'S'*+^_ can only increase with 
increasing 5*^^. 

It is instructive to understand better the results in the scalar scenario. Inspecting the 
formulae for the Wilson coefficients we arrive at an important relation: 

= SZ' - TT, (131) 

where the shift is related to the minus sign difference in the Z' and scalar propagators. 



But as seen in (84) the three observables given there, all depend on 2(ps, implying that 
from the point of view of these quantities this shift is irrelevant. As different oases 
correspond to phases shifted by vr this also explains why in the scalar case the results in 
different oases are the same. That the branching ratio can only be enhanced follows just 
from the absence of the interference with the SM contributions. In order to understand 
the signs in S** one should note the minus sign in front of sine in the corresponding 



formula. Rest follows from (131) and Table |6| 

In Fig. [g] we plot vs S^^ for A^ and cases. We observe that for Mh = 1 TeV, 
even for S*^^ significantly different from zero, A\y does not defer significantly from 
unity in both scenarios. Larger effects have been found in the Z' case as seen Fig. 4 
of 2). 

In Fig. [7] we show how the plots in Figs. |4] and [5] change when the exchanged particle 
has both scalar and pseudoscalar couplings to muons with 

Af = 2Af = 0.4. (132) 

We observe that while the correlation between S^^+^^ and S^^ is relative to A^ case prac- 
tically unmodified, the correlation between 5*^0 and B{Bs — )■ /i"^/i~) is visibly modified 
for B{Bs —7- yU^yU^) below the SM value while less if an enhancement is present. 
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Figure 7: S^^+^- versus S^^ (left) and S^^ versus B{Bs — ?■ fi^fi ) (right) for Mh 



The 



1 TeV in LHSl for two oases and the mixed and case with A 
magenta region corresponds to the U{2)^ limit for LHSl and the cyan region for LHS2. 
Gray region: exp la range B{Bs — )• jj-^jj.^) = {3.2^i^) ■ 10^^. Red point: SM central 
value. 



Clearly for a fixed Mh the results presented so far depend on the choice on muon 
couplings made by us. In Figs. |8] and [9] we show the corresponding plots when the lepton 
couplings Ap and A^ are varied independently in the range 0.02 — 0.04. Evidently, the 
allowed regions are now larger but the general pattern of correlations remains. These 
results are presented here only for illustration and we will not discuss this mixed scenario 
for other meson systems. 



8.1.2 The Bd Meson System 

We begin by searching for the allowed oases in this case. The result is shown in Fig. [TO] 
The general structure of the discrete ambiguities is as in the Bs case but now as expected 
the selected oases in SI and S2 differ significantly from each other. In fact this figure 
has the same phase structure as Fig. 6 in [l] except that the allowed values of Sis are 
reduced with respect to the Z' case for the same reason as in the Bg system: the relevant 
hadronic matrix elements are larger. 

Let us first concentrate on S2 scenario and the case. Our colour coding is such that 

• In the general case yellow and green allowed regions correspond to oases with 
small and large 613, respectively. 

• In principle in the U{2)^ symmetry case we could again show the reduced re- 
gions with magenta and cyan for LHSl and LHS2, respectively but this reduction 
amounts typically to 5 — 10% at most and it is more transparent not to show it. 
This small impact of U{2)^ symmetry in the B^ system is evident from Table ^ 
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Figure 8: Correlation plots as in Fig. and[^ but now with lepton coupling Ap G 

[0.02,0.04], A^= (left) and As G [0.02,0.04], Ap = (right) m LHSl. Gray region: 
exp la range B{Bs — > = (3.2;^|2) • 10^^. Red point: SM central value. 



In the right panel of Fig. 11 we show S^Ks B{Bd — )■ yU^/i ). This result should be 
compared with the one for Z' case shown in the right panel of Fig. 7 in [l]. 

In order to understand the differences between these two scenarios of NP we again look 
at the phase of the correction to P which now is given as follows: 

71 

rz'^TAo, Sz'=Si3-f3, Sao = Sz' - -. (133) 
Note that this time the phase of Vtd enters the analysis with /3 ^ 19° and /3 ~ 25° for 
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Figure 9: Correlation plots as in Fig. but now with Ap E [0.02,0.04] and As G 
[0.02,0.04]. Gray region: exp la range B{Bs — )■ fi~^fJ^~) = {3.2til) ■ 10~^. Red point: 
SM central value. 
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Figure 10: Ranges for AMd (red region) and S^Ks (blue region) for Mh = 1 TeV in 
LHSl (left) and LHS2 (right) satisfying the hounds in Eq. (121). 



SI and S2 scenario of \Vub\, respectively. We find then that in scenario S2 the phase 6z' 
is around 115° and 295° for yellow and green oases, respectively. Correspondingly ^^^o 
is around 25° and 205°. We summarize the ranges of 5^/ and in Table [6j 

With this insight at hand we can easily understand the plots in question noting that 
the enhancements and suppressions of B{Bd — )■ /^"""/i") are governed by the cosine of 
the phase of the correction: 

• In the case B{Bd — ?■ /^"""/i") is enhanced in the yellow oasis but suppressed in 
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Figure 11: S^Ks versus B^B^ — )■ /^^/^^) scenario for Mh = 1 TeV in LHSl (left) 

and LHS2 (right) in the yellow and green oases as discussed in the text. Red point: SM 
central value. 



the green oasis. 

• In the Z'-case the behaviour is opposite: Bi^B^ — )• /i^/i^) is suppressed in the 
yellow oasis but enhanced in the green oasis. 

• For the choice of parameters NP effects are a bit larger in the case. 

Note that in both cases the requirement on S,^Ks and AM^ forces B{Bd — )■ /i^yu^) to 
differ from the SM value. In the case these enhancements and suppressions amount 
up to ±70% for Mh = 1 TeV. They increase with decreasing S^Ks- 

Note that because of the correlation between B{Bd — )• fi^fi^) and S^Ks ^i-nd the fact that 
the latter is already well determined, the range of 613 cannot be large and consequently 
also the ranges for the phases 6z' and 5^0 are small allowing thereby the identification 
of the right oasis by B{Bd — )■ /^"^yU") only. 

While the correlation between B{Bd — )■ /i^/i^) and S^Ks offers a distinction between 
Z' and pseudoscalar scenario even more interesting in this respect is the correlation 



between B{Bd l^^fJ' ) and S'|^+^_. We show this correlation in Fig. 12 (right panel). 
Comparing this result with right panel of Fig. 8 in [l] we observe that in LHS2 these two 
observables are correlated within A'^ case but anticorrelated in the Z' scenario. Also 
this behaviour follows directly from the phase structure of NP contributions in these 
two cases. Measuring only the signs of shifts in these two observables relative to the 
SM values will uniquely tell us which of these two NP scenarios could be at work and 
which not. 

We next turn to LHSl scenario for \Vub\- We observe that the phase 613 is lower than in 



the case of scenario S2 but S13 is basically the same. Using (133) we can again calculate 



the phases of NP physics contributions to P. We find that now 6z' is around 70° and 
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Figure 12: — )■ /i/i) versus 5"^+^- case /or M^/ = 1 TeV in LHSl (left) and 

LHS2 (right) for the green and yellow oases as discussed in the text. Red point: SM 
central value. 



250° for yellow and green oases, respectively. Correspondingly is around —20° and 
160°. We summarize the ranges of 5z' and in Table [6j 



With this insight at hand we can easily understand the plots in the left panels in Figs. 11 



and 12 and analogous plots in the left panels of Figs. 7 and 8 in [l]. In particular as 
there is a flip in sign of •5'^+^- when moving from LHSl to LHS2 while 



ig.[T|t 

o quali 



seen in Fig 
there is no 
found in the Z' case. 



itative change in the case of B{Bd /U+yU ). The opposite behaviour is 



Therefore what distinguishes LHSl from LHS2 in both NP scenarios is the sign of 
the correlation between S'^+^^ and B{Bd — )■ In the A° positive 5'';+ 

implies suppression of B{Bii — )■ fi^fj^^) in LHSl but enhancement in LHS2. Note that 
this pattern is independent of the sign of Ap^ coupling as this coupling enters both 



observables. On the other hand the flip of this sign would interchange colours in Figs. 11 



and 12 As seen in Figs. 7 and 8 in p!] in the Z'-case the behaviour is opposite to the 



one found in the case: anti-correlation in LHSl and correlation in LHS2 between 
S'|^+^_ and B{Bci — > /i^/U^) in the case is changed respectively to correlation and 
anti-correlation in the Z'-case. This means again that once we will know whether LHSl 
or LHS2 is chosen by nature the measurements of S'^^+^j,- and B{Bd — ?■ /i'^yu") will tell 
us which of the two NP scenarios are favoured. However, we are aware of the fact that 
while the measurement of >S'*+^_ is extremely difficult, the measurement of >S'^+^_ will 
require heroic efforts and it may take decades to realize such a measurement. 



We next move to consider the case and show the results in this case in Figs. 13 



and 14 that correspond to Figs. 11 and 12 in the A^ case, respectively. We observe the 



following differences between A^ and cases: 

• As in the case of Bs — )■ fi'^fJ'~ there is no dependence on the oasis considered and 
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Figure 13: S^Ks versus B{Bd — )■ /i^/i^) m scenario for Mh = 1 TeV in LHSl 
(left) and LHS2 (right) in the yellow and green oases that overlap here. Red point: SM 
central value. 



in all plots in Figs. 13 and 14 -6(5^ — )■ jj,'^ ji ) is always enhanced as opposed to 
the case where it could also be suppressed. 

The asymmetry S'l^+^. is negative and positive in LHSl and LHS2, respectively, 
while in the case both signs were possible in LHSl and LHS2. 



In order to understand the signs of 5'^+^- in this case one should just use Table 6 and 
the relation (131) in the Bd system. Effectively these signs in Figs. 13 and ll can 



be obtained from the corresponding plots in Figs. 11 and 12 by simply removing the 
regions with suppression of B{Bct — )■ Therefore the distinction between S and 

P will only be easy if this branching ratio will turn out to be suppressed with respect 
to its SM value. 



8.1.3 The U{2Y Limit 

We have investigated how the results presented until now are modified when the flavour 



f/(2)^ symmetry 74 80 is imposed on the H couplings. As pointed out in jo] in this 
case ifB^ = ^PBs which in turn implies not only the correlation between CP asymmetries 
S^Ks s^nd but also a triple S^Ks ~ Si)4> ~ \ Vub\ correlation. 

Usually, when considering the case of U{2)^ broken by the minimal set of spurious, the 
Mf/(2)^ case, only SM operators are involved. Yet, as mentioned in M in the case 



of 2HDM^jpY 81 with flavour blind phases dominantly in the Higgs potential and the 
dominance of scalar left-handed currents also the U{2)^ structure of scalar contributions 
to FCNC transitions is obtained. Thus only the LHSl and LHS2 scenarios are involved 
in this case. In what follows we will confine our discussion to Bg and B^ systems as NP 
effects in the K system are much less interesting in the NP scenarios considered in this 
paper. General discussion can be found in (ol. 
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Figure 14: B{Bd — )• fifl) versus 3'^^+^- in case for Mh = 1 TeV in LHSl (left) and 
LHS2 (right) for the green and yellow oases (they overlap here) as discussed in the text. 
Red point: SM central value. 



Until now NP effects in the observables in and systems where uncorrelated but 
now they are correlated with each other due to the relations: 



_Sl3_ 



td\ 



S23 



6i3 - 6' 



23 



(134) 



ts 



Thus, once the allowed oases in the B^ system are fixed, the oases in Bg system are 
determined. Moreover, all observables in both systems are described by only one real 
positive parameter and one phase, e.g. (§23,^23)- 

We also have in this case [qI 



S 



s- 



(135) 



for which formulae in the H and A case can be found in Subsection 4.5 



In Fig. 15 we combine Figs. |2]and 10 using the ?7(2)^ symmetry relations in (134). In 
the U{2)^ limit the allowed oases get smaller. This decrease turns out to be not very 
pronounced in the case of (sis, S13) oases as they were already small as seen in Fig. 10 
but has a significant impact on (§23, 523) oases which where much larger as seen in 
Fig. |2} Moreover the fact that the results in the B^ system depend on whether LHSl or 
LHS2 is considered is now transfered through the relations in (134) into the Bs system. 



This is clearly seen in Fig. [15} in particular the final oases in cyan in LHS2 are visibly 
smaller than the magenta oases in LHSl due to the required shift of S^Ks- 

Inspecting the ranges for the phases in the last four rows of Table [6| which now apply to 
both Bg and B^ systems and comparing them with the remaining rows of these tables 
we can predict the impact of the imposition of U{2Y symmetry on the results in both 
systems presented so far: 

• In the case of Bg system in which previously there was no distinction between 
LHSl and LHS2 the changes are as follows. The plots in Figs. |4] and |5] still apply 
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Figure 15: Ranges for AMg (red region), S^^ (blue region), AM^ (green region) and 
S^Ks (yellow region)for = 1 TeV in LHSl (left) and LHS2 (right) in the U{2)^ 
limit satisfying the bounds in Eq. (126) and (121). The overlap region of all four 
regions is shown in magenta in LHSl and in cyan in LHS2. 



but the allowed regions get smaller and some of them are valid only for LHSl 
{magenta) and other for LHS2 (cyan). We also note that the imposition of U{2)^ 
symmetry favours regions away from the SM point. 



• In the Bd system the plots in Figs. [Tl] - [14] have the same structure as previously 
and as stated previously the effect of the imposition of U{2)^ symmetry is so 
small that we do not show magenta and cyan areas in this case. They would 
cover almost completely the yellow and green areas, respectively. 

These expectations could be already made by comparing the right panel of Fig. 3 and 
Fig. 25 in [l] in the case of Z' scenario. Indeed also there the imposition of U{2)^ 
symmetry reduces the allowed ranges significantly in the Bg system and make a clear 
distinction between LHSl and LHS2 scenario which was absent previously. We see again 
how important the determination of \Vub\ is. Knowing future precise values of \Vub\ as 
well as 5*^0 and B{Bs,d — ^ f^~^f^~) will confirm or rule out this scenario of NP. These 
correlations are particular examples of the correlations in MU{2)^ models pointed out 
in [9]. What is new here is that in a specific model considered by us the \Vub\ — S^^ 
correlation has now also implications for B{Bs — > ^'^ ^jT) and S"*^. 



In Fig. 16 we show B{Bci — )■ versus B{Bs — ?■ ^^fi") for the ^4° and cases. As 

expected on the basis of a general discussion in [o] there is a very strong correlation 
between these two branching ratios. Again, while in the case both branching ratios 



8 An Excursion through H Scenarios 



49 



LHS & 4" 



LHS & H*' 




Figure 16: B{Bd — )■ fi~^fi ) versus B{Bs — )■ /i"*"// ) /or Mh = 1 TeV in case (left) 
and case (right) in the U{2)'^ limit satisfying the bounds in Eq. (126) and (121). 



can be enhanced or suppressed with respect to the SM, they can be only enhanced in 
the case. 



8.2 The RHSl and RHS2 Scenarios 
8.2.1 First Observations 

We will now investigate H scenario with exclusively RH couplings to quarks. Now in the 
RHSl and RHS2 scenarios only RH couplings to quarks are present in H contributions. 
As QCD is parity conserving, the hadronic matrix elements for operators with RH 
currents as well as QCD corrections remain unchanged. The expressions for AF = 2 
observables in RHSl and RHS2 scenarios as well as the corresponding constraints have 
precisely the same structure as in the LHSl and LHS2 cases just discussed. Therefore 
the oases in the space of parameters related to RH currents are precisely the same 
as in LHSl and LHS2 scenarios, except that the parameters s^- and 6ij parametrize 
now RH and not LH currents. Anticipating this result we have not introduced separate 
description of LH and RH oases. Yet, in the case of AF = 1 observables several changes 
are present which allow in principle to distinguish the RHSl and RHS2 scenarios from 
the corresponding LHSl and LHS2 scenarios. 



8.2.2 The Bs Meson System 

On the left in Figs. |4]and jsjwe have shown 5'^+^- vs S^^ in the LHSl scenario for 
and cases, respectively. Analogous plots for the correlation of S^^ vs B{Bs — )■ 
in the LHSl scenario are shown in Figs. |4] and [5j Inspecting the related formulae for 
RHS scenario we conclude similarly to the Z' case that in the case of A° these plots are 
also valid for RHSl scenario except that the colours should be interchanged. Therefore 
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on the basis of observables considered here it is not possible to distinguish between 
LHSl and RHSl scenarios because in the RHSl scenario one can simply interchange 
the two oases to obtain the same physical results as in LHSl scenario. We also note 
that reversing simultaneously the sign of Ap would keep also the oases unchanged. 

The situation is even simpler in the case of case. As the plots in question did not 
depend on oasis considered, the correlations in this case are identical in RHSl and LHSl 
independently of oasis considered. 

Clearly as in the LHSl scenario this result represents a test of the RHSl scenario but 
if one day we will have precise measurements of S^+^^, S^^ and B{Bs — )■ /i^/U^) we will 
still not be able to distinguish for instance whether we deal with LHSl scenario in the 
blue oasis or RHSl scenario in purple oasis. 

In principle, one could make a distinction between LHS and RHS scenarios by consid- 
ering model independent bounds from B — )• and B — )■ K*fi~^fi~ on the Wilson 



coefficients of the scalar operators. However, as discussed in Subsection 8.5, this is 



presently not the case. This should be contrasted with Z' analysis in yj where in fact 
such a distinction could be made. 



8.2.3 The Ba Meson System 

Similarly to the Bg case the structure of oases is as in Fig. 



10 Moreover, the results 



in in Figs. 11 -14 are valid for RHSl and RHS2 scenarios except that in the A case 



the colours should be interchanged, while there is no modification in the case. 
Thus we cannot distinguish between LHS and RHS scenarios on the basis of considered 
observables. Clearly the study of 6 — )■ dfi'^fi~ transitions could help in this context but 
they are more challenging both theoretically and experimentally. 



8.3 The LRSl and LRS2 Scenarios 
8.3.1 First Observations 

If both LH and RH currents are present in NP contributions, the pattern of flavour 
violation can differ from the scenarios considered until now in a profound manner. If 
the LH and RH couplings differ from each other, the number of parameters increases and 
it is harder to get clear cut conclusions without some underlying fundamental theory. 
On the other hand some of the "symmetries" between LHS and RHS scenarios identifled 
above are broken and the effect of RH currents in certain cases could in principle be 
better visible. 

Here in order to keep the same number of parameters as in previous scenarios we will 
assume a left-right symmetry in the if-couplings to quarks. That is the LH couplings 
Al are equal in magnitudes and phases to the corresponding RH couplings A/?. In this 
manner we can also keep the same parametrization of couplings as in previous scenarios. 
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Before entering the details let us emphasize two new features relative to the cases in 
which either LH or RH couplings in NP contributions were present: 

• NP contributions to AF = 2 observables receive now new LR operators, whose 
contributions are enhanced through renormalization group effects relative to SM 
operators, however as scalar LL and RR operators are also enhanced by such 
effects the difference between LL (RR) scenarios and LR scenario in the scalar 
case is much smaller than in the Z' case. 

• NP contributions to 5^ ^ — )■ fi'^fi" and Kl — )■ fi'^fi~ vanish eliminating in this 
manner 5*+ _ and B{Bs^d /^^A* ) a-s basic observables in the identification of 
acceptable oases. On the other hand B — > K*fi^^~ and B — )■ Kfi~^fi~ receive still 
NP contributions and can help in this context. 

While >S'*+^_ cannot help in the identification of the optimal oasis in the LR scenarios 
they are non-vanishing: 

Sl^^^ = -sm{2^B,). (136) 
While rather small they offer a clean test of the LR scenarios. 



8.3.2 The Meson System 

We begin the search for the oases with the Bg system proceeding with input parameters 
as in the previous scenarios. The result of this search for Mh = 1 TeV is shown in 



Fig. 17, where we show the allowed ranges for (s23, ^23). The red regions correspond to 
the allowed ranges for AMs, while the blue ones to the corresponding ranges for S^^. 
The overlap between red and blue regions identifies the oases we were looking for. 

The notations are as in previous cases but it should be kept in mind that the parameters 
(§235 ^23) describe both LH and RH couplings. 



In order to understand the structure of oases in Fig. 17, that differs from the ones found 
so far, we note that the matrix element of the dominant (^2^ operator has the sign 
opposite to the dominant Qf^^ operator. Therefore, this operator naturally suppresses 
AMg with the phase ^23 centered in the ballpark of 0° and 180°, that is shifted down 
by roughly 90° relatively to the LHS scenarios. As the matrix element of (^2^ is larger 
than that of Qf^^ operator in LHS and RHS scenarios, §23 has to be sufficiently smaller 
to agree with data. 

The crucial role in the Bg meson system in this scenario, in the absence of NP contri- 
butions to Bg^d — > decays, is now played by i? — )■ K*fi~^fi' and B — t- Kfi^^~. We 
will discuss the latter decays at the end of this section. 
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AM, & S.,,/, LRS1 




0.000 0.001 0.002 0.003 0.004 0.005 0.006 0.007 
S23 



Figure 17: Ranges for AMg (red region) and S^^ (blue region) for Mfj 
LRSl satisfying the hounds in Eq. (126). 



1 TeV in 



8.3.3 The Ba Meson System 



The structure of oases in this case is given in Fig. 18, As we do not have B{Bd — )■ /^^/^~) 
to our disposal and h — j- di^i~ decays are challenging this system is not very useful to 
provide tests of LRS scenarios without some fundamental theory. 

Due to the sign of the matrix element of the dominant (^2^ operator in both LRSl and 
LRS2 the mass difference AM^ is naturally suppressed. The requested size of this sup- 
pression together with significant suppression of S^^Xg in LRS2 and slight enhancement 
of it in LRSl governs the structure of the phases. 



8.4 The ALRSl and ALRS2 Scenarios 



We include this case as well because it has not been discussed in the literature but it is 
an interesting NP scenario for the following reasons: 

• NP contributions to AF = 2 observables are dominated as in LRS scenarios by 
new LR operators but as the sign of LR interference is flipped some differences 
arise. 

• NP contributions to Bd^s A*^/^^ enter again with full power. Therefore these 
decays together with 5'^+^- offer as in the LHS and RHS scenarios some help in 
the identification of acceptable oases and to study differences between scalars, 
pseudoscalars and Z' bosons. 
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AMrf a S^„, LRS1 



AMtf S S^Kg. LRS2 
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Figure 18: Ranges for AM^ (red region) and S^Ks (blue region) for Mh = 1 TeV in 
LRSl (left) and LRS2 (right) satisfying the hounds in Eq. (127). 



• The phase structure of the oases is as in LHS scenario but due to enhanced 
hadronic matrix elements of LR operators the mixing parameters Sij are decreased. 

• NP contributions to — )■ vr^f+f" vanish in this scenario. 

In view of this simple structure of modifications with respect to LHS scenario, all plots 
have the same structure as LH scenarios but NP effects are smaller. Therefore we will 
not show these plots. 



8.5 Implications of 6 ^ si^i Constraints 

Presently the NP effects found by us are consistent with the experimental data on 
Bs,d yU+Z^"- However, also the data on B XJ+i-, B K*i+i- and B Ki+r 
recently improved a lot and it is of interest to see whether this has an impact on 
our results. It should be emphasized that B — t- Xgi'^i", B — )■ K*i~^i~ are not as 
theoretically clean as Bg — )■ fi^fi^ because of the presence of form factors. However in 
the case of i? — i- K£^£^ progress in lattice calculations of the relevant form factors is 
expected soon and as stressed in particular in fT6^ a simultaneous consideration of this 
decay together with Bs — )■ provides useful tests of extensions of the SM. Indeed, 

while Bg — )■ is sensitive only to the differences Cp — C'p and Cs — C^, the decay 

B — > Ki~^i~ is sensitive to their sums Cp + Cp and Cs + C'g. A very extensive model 
independent analysis of Cp{C'p) and Cs{C'g) in the context of the data on Bg — t- fJ^'^fi" 
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and B — t- Kl^l has been performed in 16 . Adjusting their normahzation of Wilson 

(137) 



coefficients to ours the final result of this paper reads: 



which implies 



mfelC^'^l < 0.7, mb|C);^|<1.0 
1 1 AtJH)^f{H) 



1 1 Ai\iH)AfiH) 



V*Vt 



ts^tb 



or equivalently for Mh = 1 TeV 

\~S2^Af{H)\ < 0.00115, 
|523Af < 0.00164. 



< 0.70 



< 1.0 



(138) 
(139) 

(140) 
(141) 



The largest values of S23 used in our analysis are 0.0041 in LHS and RHS scenarios 
with smaller values for LR and ALR scenarios. We find then 0.000164 and 0.000082 for 
the two products respectively. Therefore, these bounds do not have any impact on our 
results. 



9 The K Meson System and SM Higgs 
9.1 The K Meson System 

As we already stated previously we do not expect any visible effects in — )■ -k^vv and 
Kl — )• 7r°z/z/ and our discussion will concentrate on sk, Kl — )• fi~^fi~ and Kl — )■ ir^i'^i' . 



As seen in ( 128 ) the constraints from AF = 2 observables are weaker than in previous 



cases. Yet as seen in Fig. 19, obtained within LHSl and LHS2 scenarios, it is possible 
to identify the allowed oases. These plots have the same phase structure as the plot in 
Fig. 9 of [T] for Z' scenario except that S12 is by a factor of five smaller because of the 
enhanced matrix element of the relevant scalar operator. 

Due to weaker constraints in the K system the oases are rather large. We have two 
oases in SI: 

Ci{Sl) : 0° < 612 < 90°, C2{S1) : 180° < 612 < 270° (142) 
and only one oasis in S2: 

Ci{S2): 0°<5i2< 360°. (143) 

With these constraints at hand we have calculated the branching ratios for — ?■ /J'^IJ'" 
and Kl — )■ TT^i^i~ decays. We can summarize our results as follows: 
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AM If & 1*^1, LHS1 



AWI« & le^l, LHS2 




3 




^ 3 



0.00002 0.00004 O.OOO06 0.00008 0.00001 0.00002 0.00003 0.00004 0.00005 

S12 Si2 

Figure 19: Ranges for AM^ (red region) and Ek (blue region) (LHSl: left, LHS2: 



right) for Mh = 1 TeV satisfying the bounds in Eq. (128) 



1.40 



LHS1 & A° 



_ t.6r 
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E{Ki_^ti* (i-] [10-^1 



1.2^ 
0.0 



LHS2 & A° 




0.5 1.0 1.5 



2.0 



Figure 20: B{Kl — > 7r°/i+/i ) versus B{Ki — ?■ ) for P scenario and LHSl (left), 
LHS2 (right). Red point: SM central value. 



• NP effects in — )■ n^e^e are totally negligible both for and cases. 

• NP effects in — n^fi^fi' are larger but amount to at most ±5% at the level 
of the branching ratio which is also negligibly small in view of large theoretical 
uncertainties. 

• The short distance branching ratio for — t- fi^fi can a be modified up to ±50% 
in the case. Still such effects are fully consistent with the upper bound on this 
branching ratio. NP effects in the H'^ case are much smaller. 
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These results are rather disappointing but allow to distinguish scenarios discussed here 
from Z' scenario, where effects have been found to be larger. As an example we show 



in Fig. 20 the correlation between B{Kl — )■ n^fi'^fi ) and B{Kl — )■ /i"*"/! ) in LHSl and 
LHS2. 

The effects in other scenarios are rather uninteresting as well and we will not present 
any results for rare K decays in them. 

However, it is of interest to see how the oases change in the presence of LR operators. 



This we show in Fig. 21 Due to the presence of LR operators the structure of oases is 
different than in LHSl and LHS2 scenarios. While the shape of the single oasis in the 
LRS2 case is similar to the LHS2, for LRSl the oases are shifted by 90°: 

Ci{Sl) : 90° < 6i2 < 180°, C2{S1) : 270° < 612 < 360°. (144) 

NP effects in Kl — )■ fJ'~^P-~ vanish in these scenarios and in Kl — )■ ir^i'^i^ they are 
negligible. Therefore we do not show any plots. 

Concerning NP contributions to K ^ rruu decays all scalar scenarios could turn out 
one day to be interesting if the data on observables in Bg and Bd systems will show the 
presence of NP but negligible NP effects in K ^ ttuu. 



9.2 Flavour Violating SM H Boson 



We will next turn our attention to flavour violating couplings of the SM Higgs (h) that 
can be generated in the presence of other scalar particles and or new heavy vectorial 
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fermions with +2/3 and —1/3 electric charges. In the case considered by us, new 
quarks with —1/3 charges are essential for generating flavour violating couplings to SM 
down-quarks but the presence of heavy quarks with +2/3 charges could be relevant for 
charm physics. Moreover, such heavy fermions could contribute to rare K and B decays 
through loop diagrams. In what follows we will not consider these loop contributions 
as they would lead us beyond the scope of our paper. 

The strategy and formalism developed in the previous sections can be used in a straight- 
forward manner for the case of h flavour- violating couplings to quarks. Using the general 
relations of Section |6] one then finds that these couplings have to be significantly smaller, 
by roughly an order of magnitude relatively to the corresponding couplings of a heavy 
scalar with a mass in the ballpark of 1 TeV. Still the presence of such contributions 
can remove all possible tensions within the SM in AF = 2 transitions without being in 
conflict with constraints from rare decays, where the SM Higgs contributions, in spite 
of a low Higgs mass, are negligibly small. 

The reason for the smallness of SM Higgs contributions to rare decays originates from 
the smallness of scalar Higgs coupling to with A^^(/i) = 1.2 x 10~^. It is roughly a 

factor of 40 smaller than the largest scalar coupling allowed by for a scalar 

with Mh = 1 TeV. As the correlations between AF = 2 and AF = 1 transitions in 
Section [6] show, this smallness of the muon couplings of SM Higgs can be compensated 
partly by the smallness of its mass. But this compensation is insufficient to make the SM 
Higgs contributions to rare B^^d and K decays of any relevance as these contributions 
do not interfere with the SM contribution from Z-penguins and box diagrams. They 
are suppressed by the square of Ag^{h). 

10 Summary and Conclusions 

In this paper we exhibited the pattern of flavour violation in models in which NP effects 
are dominated by tree-level heavy pseudoscalar (A) or scalar (H) exchanges under the 
assumption that the theoretical and experimental errors on various input parameters 
will decrease with time. In particular we have identified a number of correlations 
between AF = 2 and AF = 1 processes that will enable in due time to test this NP 
scenario. Our detailed analysis of these correlations in Section |8] shows that in the Bs 
and Bd systems a very rich pattern of NP effects is present while this is not the case 
in K decays. This is partly opposite to Z' and Z scenarios considered in [l] where the 
largest effects have been found in the K system, even for masses of Mz' outside the 
LHC reach. 

Our results are summarized in a number of plots that have been obtained in various 
scenarios for the H couplings and for inclusive and exclusive values of \Vub\- We list 
here only few highlights: 

• For each scenario we have identified allowed oases in the parameter space of 
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the model. In each oasis particular structure of correlations between various 
observables will in the future either favour or exclude a given oasis. 

For the near future the correlations involving S^Ks^ and B{Bd — ?■ yu^/i^) and 
B{Bs — 7- fi^fJ'~) will be the most interesting as the data on these four observables 
will be improved in the coming years, sharpening the outcome of our analysis and 
possibly ruling out some oases and scenarios of the couplings. 

Most importantly we have found that various correlations involving 5*^^, B{Bs — t- 
f^~^fj'~) and the CP asymmetry 5'^+^- show profound differences between the sce- 
nario with a pseudoscalar tree-level exchange and Z' exchange. As we explained 
in detail, these differences are directly related to the difference in the fundamental 
properties of the particles involved: their spin and CP-parity. As far as the last 
property is concerned also differences between the implications of the pseudoscalar 
and scalar exchanges have been identified. In particular the scalar contributions 
can only enhance B{Bs — )■ and are invariant under the interchange of 

two oases in parameter space involved, which is not the case of pseudoscalar ex- 
changes where the branching ratio can also be suppressed. The symphony of 



plots in Subsection |8.1[ where the most important results of our paper are shown 
will be helpful in monitoring further developments in the measurements of the 
observables in question. 

Analogous comments apply to the correlations involving S^Ksj B{Bd — )■ 

and the CP asymmetry 5*^+^- for which the symphony of plots, also presented in 

Subsection 8.1, is doubled in view of the dependence on the value of \Vub\- 



We have demonstrated that the imposition of f/(2)^ symmetry on pseudoscalar 
and scalar quark couplings has a profound impact on the correlations in the Bs 
system with much smaller effects in the B^ system. 

We have also pointed out additional differences between Z' and pseudoscalar or 
scalar tree-level contributions related to channels with neutrinos in the final state, 
where in the Z' case these contributions could be very large but are expected to 
be negligible in NP scenarios considered here. 

Our short analysis of flavour-violating SM Higgs-couplings shows that in the case 
of rare Bg^d and K decays, the SM Higgs contributions are irrelevant due to the 
smallness of the Higgs coupling to muons. On the other hand such contributions 
could in principle remove all tensions within AF = 2 observables observed within 
the SM. 



We close our paper by Fig. 22 in which we show the correlations involving 
and B{Bs — )■ Z^"^/^") combining information of Fig. [sjfor the tree-level scalar and pseu- 
doscalar exchange and include also tree- level Z' exchange. The lepton couplings are 
not fixed but varied in the following ranges: 



e [0.02,0.04] and A7(Z') 
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Figure 22: Overlay of the correlations for S^^ versus S^^ (top left), versus S^p^ (top 
right) and S^^ versus B{Bs /U^/U^) (bottom) for tree level scalar (cyan), pseudoscalar 
(red) and Z' (blue) exchange (both oases in same colour respectively) in LHS. The lepton 
couplings are varied in the ranges \/S.^g'p{H)\ E [0.02,0.04] and A'^{Z') G [0.3,0.7]. 



[0.3,0.7]. Further we do not distinguish between the two different oases here. The 
patterns aheady identified previously and summarized above are clearly visible in these 
plots. 

We are aware of the fact that some of the correlations presented by us would be washed 
out if we included all existing uncertainties. Yet, our simplified numerical analysis 
had as the main goal to illustrate how the decrease of theoretical, parametric and 
experimental uncertainties in the coming years might allow to exhibit certain features 
of NP, even if deviations from the SM will be only moderate. In this manner we 
have uncovered a world of correlations present in NP scenarios, where new effects are 
dominated by flavour- violating couplings of a heavy neutral pseudoscalar and scalar. 
In fact, within the coming years the size of the assumed uncertainties in our analysis 
could likely become reality not only because of improved experimental data but also 
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improvements in theory, in particular lattice calculations of hadronic matrix elements, 
Bi parameters, form factors and weak decay constants. 

We are looking forward to improved experimental data and improved lattice calcula- 
tions. The correlations identified in this paper will allow to monitor how simple NP 
scenarios discussed by us face the future precision flavour data. 
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